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ON THE USE OF THE INDICIAL FUNCTION CONCEPT IN THE ANALYSIS OF UNSTEADY 

MOTIONS OF WINGS AND WING-TAIL COMBINATIONS 

By Mtjkkat Tobak 


SUMMARY 

The concept of indicial aerodynamic functions is applied to 
the analysis of the short-period pitching mode of aircraft. By 
the use of simple physical relationships associated with the 
indicial-function concept, qualitative studies are made of the 
separate effects on the damping in pitch of changes in Mach 
number, aspect ratio, plan-form shape, and frequency. The 
concept is further shown to he of value in depicting physically 
the induced effects on a tail surface which follows in the wake of 
a starting forward surface. Considerable effort is devoted to 
the development of theoretical techniques whereby the transient 
response in lift at the tail to the wing wake may he estimated. 
Numerical results for several representative cases are presented, 
and these are analyzed to reassess the importance of the contri- 
bution to the rotary damping moment of the interference lift at 
the tail. 

INTRODUCTION 

In tlie classical study of the longitudinal motion of an 
aircraft, it is usually found that the motion resulting from a 
small equilibrium-destroying disturbance consists of two 
modes: one, a lightly damped, low-frequency motion at 
essential^ constant attitude, called the phugoid oscillation; 
the other, a rotary-pitching and plunging oscillation of high 
frequency (relative to the phugoid frequency) called tlie 
short-pci-iod oscillation. The phugoid oscillation has gen- 
erally been described as resulting from a slow interchange of 
potential and kinetic energy as the aircraft experiences 
periodic variations in airspeed and altitude. The character 
of the phugoid .motion as influenced by airspeed, altitude, and 
aircraft geometry has been well understood for some time 
(see, e. g., ref. 1). The short-period motion, on the other 
hand, having in the past been found to be highly damped and 
of short duration, has been tUe cause of no concern. Its 
characteristics therefore have not been as fully investigated 
as those of the phugoid oscillation. As flight speeds have 
progressed to the transonic and supersonic domains, however, 
the loss of rotary damping occurring in practically all aircraft 
at speeds near the sonic speed has caused renewed interest 
in the short-period pitching mode. Unlike the easily con- 
trolled phugoid oscillation, the deterioration of damping in 
the short-period mode is of serious concern to the pilot, since 
the period of oscillation can be of the same order of magnitude 
as the pilot’s reaction time. The oscillation may therefore 
be difficult or even impossible for the pilot to control man- 
ually. Furthermore, the additional load imposed upon the 


airframe due to a rapid growth of the amplitude of a nega- 
tively damped oscillation makes possible the occmrence of 
structural failure. It is therefore of considerable interest to 
obtain an understanding of the nature of the short-period 
mode, parallel to that which has been gained of the phugoid 
mode. 

One means of viewing the aerodynamic phenomena occur- 
ring during the short-period oscillation from a fundamental 
standpoint is through application of the concept of indicial 
functions. In this approach, the variations with time of the 
airci'aft angle of attack and angular velocity are replaced by 
a large number of small instantaneous or step changes. The 
transient aerodynamic reactions to these step changes are 
termed “indicial functions” and have been calculated for 
several classes of isolated wings (refs. 2 to 6). By suitable 
superposition of these results (refs. 7 to 9) the aerodynamic 
forces and moments caused by the given maneuver can be 
studied. It will be the primary purpose of this report to 
make such a study for the si.mplified case of an aircraft per- 
forming single-degree-of-freedom rotary oscillations, a ma- 
neuver which corresponds to the short-period oscillation when 
the plunging velocity of the aircraft is zero. To effect this 
end systematically the report is organized in two main cate- 
gories: First the motion of a tailless aircraft is studied, and 
here existing theoretical information and the use of simple 
physical relationships associated with the indicial-function 
concept enables qualitative studies to be made of the sepa- 
rate effects on the aerodynamic forces and moments of 
changes in Mach number, aspect ratio, plan form, and fre- 
quency. In the second part, the motion of a tailled aircraft 
is studied, again by use of the indicial-function concept. 
Here, however, additional theoretical information is required 
to account for the interference effect of the wing and its 
wake on the transient lift at the tad, and several sections are 
devoted to the development of theoretical techniques where- 
by this effect may be estimated. Results of computations 
based on these techniques are analyzed, and the importance 
of the contribution of the interference lift at the tail to the 
rotary damping moment is re-established. In all of the 
above work, the aerodynamic forces and moments are those 
derived from analyses of the linearized equations of potential 
flow. Thus, the usual limitations imposed by the linear 
theory on the applicability of the results are in force here 
- as well. However, in a somewhat different vein of inquiry, 
a final section is devoted to consideration of a problem in 
aircraft dynamics involving forces of nonlinear character. 
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The effect on the rotary damping moment of an aerodynamic 
restoring moment containing a common type of nonlinearity 
is studied, and conditions are pointed out under which the 
development of a self-sustained rotary oscillation may be 
possible. 

NOTATION 


A 

a. 
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. P 

aspect ratio, 


lift coefficient. 


m 

g,oS 


pitching-moment coefficient. 


pitching moment 
q„Sc 


moment of inertia 
imaginary part 

free-stream Mach number, — 

' ft 


real part 

wing (or tail) area 
flight speed 

speed of sound in free stream 
wing (or tail) span 
local chord 

wing (or tail) root chord 

forward wing mean aerodynamic chord. 


r 


(local chord )^ dy 


two-dimensional lift coefficient, 


base of natural logarithms 


qoc 


reduced frequency parameter, 


distance between moving vortex and leading edge 
of tail 

slope of wing (or tail) leading edge 
local loading (pressure on lower surface minus 
pressure on upper surface) 
angular velocity due to pitching 

free-stream dynamic pressure, ^ o 


aot' 

time 

time required following an instantaneous change 
in angle of attack or angular velocity for the 
transient lift or moment to attain steady state 


w perturbation normal velocity in plane of wing or 

tail 

Wa averaged downwash in the plane of the tail due to 

steady angle of attack of forward wing 
Wf, averaged downwash in the plane of the tail due to 

steady pitching velocity of forward wing 
z,y,z Cartesian coordinates 

Xa.c. distance from leading edge of M.A.C. to aero- 
dynamic center 

Xa distance from leading edge of M.A.C. to axis of 

rotation 

Aa;„ Xa.c.—Xo 

a angle of attack of wing center line with respect to 

free-stream direction (fig. 1) 

P VIM/-1I 

<p distance traveled, measured in half M.A.C. 

lengths, subsequent to an instantaneous change 

2V t' 

in angle of attack or angular velocity, -° 

$ perturbation velocity potential 

A<i> jump in perturbation velocity potential in plane 

of wing or tail 

d angle of wing center line with respect to horizontal 

axis (fig. 1) 

V acute angle between wing plane of symmetry and 

trailing edge (fig. 15) 
p„ free-stream density 

« angular frequency of oscillation 

(Ta distance traveled, measured in half M.A.C. 

lengths, in the time interval t', 2V„t'Jc 

When a, a, and q are used as subscripts, a nondimensional 
derivative is indicated, and this derivative is evaluated as 
the independent variable {a, a, or q) approaches zero. For 
e.xample. 




PART I— ISOLATED WINGS 

APPLICATION OF INDICIAL FUNCTIONS TO THE AERODYNAMIC THEORY OF UNSTEADY FLOWS 


One of the most useful tools in the study of unsteady 
flows is the concept of indicial aerodynamic functions, which 
may be defined briefly as the aerodynamic response of the 
airfoil as a function of time to an instantaneous change in 
one of the conditions determining the aerodynamic properties 
of the airfoil in a steady flow. Theoretical aerodynamic 
indicial functions were first derived by Wagner (ref. 2) for 
the two-dimensional wing in incompressible flow. More 
recently, these results have been extended to cover the com- 
pressible case for both subsonic and supersonic speeds (refs. 
3 and 4). In addition, theoretical indicial functions have 
now been obtained for both wide and slender triangular 


wings and rectangular wings, all for supersonic speeds (refs. 
4 to 6). 

The indicial function derives its usefulness primarily 
through the ease with which it lends itself to the powerful 
and well-established methods of the operational calculus 
(refs. 7 to 9). With the use of these methods, the aerody- 
namic forces and moments caused by arbitrary motions of 
the airframe can be studied from a fundamental standpoint. 
Because of the wide range of applicability of this means of 
approach in unsteady flow analyses, a considerable portion 
of the succeeding discussion is devoted to the fundamentals 
involved. 
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DEKNITION OF COORDINATE SYSTEM 

In the succeeding analysis the stability system of axes is 
used. The origin of the coordinate system is placed in the 
airfoil so that the y axis which is perpendicular to the vertical 
plane of symmetry is coincident with the axis of rotation of 
the airfoil; the positive branch of the x axis is pointed in the 
direction of flight; and the 2 axis lies in the vertical plane of 
symmetry, positive downward. - The angle of attack a is 
measured as the angle between the chord plane of the airfoil 
and the xy plane, and is shown as positive in figure 1. The 




angle of pitch 6 is the angle between the chord plane of the 
airfoil and the horizontal plane (an arbitrary reference) and 
is also shown positive in figure 1. Forces are measured as 
positive upward, whereas pitching moments are positive when 
tending to increase the angle of pitch in the positive direc- 
tion. When the airspeed V„ is constant, which corresponds 
to the condition under study, the translatory and angular 
motions of the airfoil with respect to any system of coordi- 
nates are defined if the time histories of the angle of attack 
a and the angle of pitch 6 and their derivatives are known. 
For purposes of clarity, two different harmonic motions of the 
aircraft are shown in figure 2, illustrating the difference 
between a flight path which involves a constant angle of 
attack and a varying angle of pitch and one which involves 
a constant angle of pitch and a varying angle of attack. 

Now consider the case of a wing executing harmonic rotary 
oscillations about the y axis while the origin of the coordi- 
nate system traverses a level path at constant velocity V„. 
This case corresponds to that of a wind-tunnel model 
mounted to permit single-degree-of-freedom rotary oscilla- 
tions, or to the short-period mode of an aircraft in flight 
when the plunging velocity of the center of gravity is zero. 


Here a and d are equal, so that the maneuver is defined by 
one variable, the time history of either a or 6. Let the angle 
of attack be a and the angular velocity be q {q=ddldt' 
=daldt'). At any instant, the normal velocity at any 
point on the airfoil surface is composed of two parts, one 
due to the instantaneous angle of attack aVo, the other due 
to the angular velocity at the same instant —qx (see fig. 3). 
These are two of the instantaneous boundary conditions of 
the unsteady flow. 

Solutions for the aerodynamic forces and moments which 
correspond to these boundary conditions may be derived by 
a number of methods involving various degrees of approxi- 
mation. In succeeding sections, the use of the concept of 
indicial functions and the principle of superposition for this 
purpose will be illustrated and compared with other current 
widely used methods. 

CONCEPT OF INDICIAL FUNCTIONS 

In order to illustrate this concept, assume that the airfoil 
under consideration has been flying a level path at zero angle 
of attack. At some time, which is designated time zero, the 
wing is caused to attain simultaneously a constant angle of 




Figure 3. — Unsteady flow boundary conditions at airfoil surface. 


attack a and angular velocity q. The normal velocity of the 
flow next to the surface of the airfoil therefore changes dis- 
continuously from zero to a pattern that is constant with 
time and identical in shape to the pattern shown previously 
in figure 3. The lift and pitching moment that result are of 
a transient character and attain their steady-state values 
corresponding to these new boundary conditions only after a 
significant interval of time has passed. It should be noted 
there exists an essential difference between the length of this 
time interval at subsonic and supersonic speeds. At super- 
sonic speeds, the vorticity shed into the airfoil wake cannot 
influence the flow about the airfoil but at subsonic speeds this 
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influence exists for nil time. The result is that the lift and 
moment reach steady-state values in a finite time at super- 
sonic speeds but approach these values asymptotically at 
subsonic speeds. In either case, however, the time responses 
in lift and moment to the step changes in a and g are termed 
“indicial functions.” Figure 4 illustrates typical subsonic 
and supersonic indicial lift responses to a step change in 
the angle of attack. 




Figure 4. — Typical indicial lift responses to step changes in angle of 

attack. 

It is obvious that the time historj'' of the wing motion 
during a short-period oscillation may be broken down into 
an infinite number of infinitesimally small step changes in 
the angle of attack and step changes in the angular velocity. 
The summation of the indicial lift and moment for these 
steps then 3 delds the total lift and moment at any prescribed 
time. In figure 5, the mechanics of the procedure are 
illustrated for an arbitrarj' angle-of-attack variation. Here, 
the given angle-of-attack variation is replaced b 3 r a number 
of small step changes. Within each step the corresponding 
response in lift is shown plotted for convenience. It is then 
apparent that the total lift at time t' is equal to the sum of 
the increments of lift in each step at time t'. As indicated 
by the leaders, however, it is clear that the increments of 



lift for the various steps at time t' are equivalent to incre- 
ments in the first step at time t' — 1[. Alternatively, then, the 
total lift at time t' can be written as 


i' A/v 

After a transformation of variables, i'— f,'=r and letting the 
increment of time approach zero, equation (1) can be re- 
written in a form of Duhamel’s integral (see, e. g., ref. 9) 

C^{t') = ^,^lCr,Sr)a{t'-r)dr (2) 


A similar procedure is carried out for the angular velocity 
variation, whereupon the total lift coefficient at the prescribed 
time t' becomes 



It should be pointed out that in this form equation (3) is 
applicable to the analysis of arbitrary motions, the only 
restriction being that the flight speed is constant. In the 
following sections, however, the application of equation (3) 
is restricted to harmonic motions having a single degree of 
freedom. The reasons for this restriction are two-fold: 
first, the motions of a statically stable aircraft in response 
to a disturbance are most generally of a harmonic nature; 
and second, such a restriction permits an assessment of the 
influence of the time rate of the airfoil motions on the 
aerod 3 mamic forces and moments. 

APPLICATION OF INDICIAL FUNCTIONS TO HARMONIC PITCHING 
OSCILLATIONS 


Consider first a pure sinusoidal pitching oscillation, the 
angle of attack being zero throughout the motion. The 
flight path for such a motion has been illustrated in figure 
2 (a). In this case, the angle of pitch is given by 

where is the maximum amplitude of oscillation and co is 
the angular frequenc 3 ''. The angular velocity is, of course, 
q{t') = d=io}9oe‘'“‘' Inserting the value for q(t') 

in equation (3) and performing the indicated operations, 
there results 

= (4) 


Note in figure 6 that (r) is equal to — F 2 (t) , 

and that F 2 (t) approaches zero as r approaches t'. 

Replacing in equation (4) b 3 =' this equality, 


Cl (t') i coc 
~W)^Wo 



(5) 


For subsonic speeds, let t' approach infinity. With this 
substitution, equation (5) thereby represents the lift coeffi- 
cient due to the harmonic pitching motion after the transient 
loading subsequent to the start of the motion has reached 
a steady periodic variation. Then separating equation (5) 
into components in phase (real part) and out of phase 


Figure 5. — Illustration of superposition process. 
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(imaginary part) with 8, there is obtained 


low frequencies, the only force of consequence during the 
pitching oscillation is the first order in frequency out-of-phase 
lift force, ikdCL^{<x^} The phase relationships for the 
harmonic pitching oscillation are indicated in figure 7. It is 
evident that the total lift leads the angle of pitch by nearly 
90°. 

APPLICATION OF INDICIAL FUNCTIONS TO HARMONIC PLUNGING 
OSCILLATIONS 

Next, consider a purely sinusoidal variation of the angle of 
attack, the angle of pitch being zero throughout the motion. 
The flight path for this motion has been illustrated in figure 
2 (b). Here, a equals where, as previously, a„ and 

o> are the maximum amplitude and angular frequency, 
respectively. Applying equation (-3) again, 

i JkB 




sin wtcIt 


( 6 ) 


Introduce the nondimensional parameters, 

r number of half M. A. C. lengths traveled in time r 


A:= 2 ^ reduced frequency. 


In terms of these parameters, equation (6) becomes, for M„ 

< 1 , 

R.P. F'2{v) cos kipdip 


LP.(^^=k F2 {(p) sin k^d<^ 


h (7) 


At supersonic speeds, equations (7) may be simplified 
somewhat since the build-up in lift is completed in a finite 
number of half M. A. C. lengths of travel Oa- In equations 
(7), therefore, the upper limits of the integrals may be re- 
placed by <7a, since beyond that point F 2 {<p) is identically zero. 


f?.P. F2{<p) cos k(fid<p 

M „>1 

I.P. (^^==k F 2 {<p) sin 


( 8 ) 


Thus it appears from equations (7) and (8) that there are 
both in-phase and out-of-phase lift forces associated with the 
harmonic pitching oscillation. Notice, however, in equations 
(8) that if the trigonometric terms are expanded and the 
reduced frequency is required to be small compared to unity 
(corresponding to the frequencies encountered in dynamic 
stability work), terms containing second and higher powers 
of k will be small compared to first-order terms. Thus, for 


ikBCi^ (ob) 


-CJ total) 


/ \ Oj, 

-ik^Bj f^{4‘)s\nk4>d4> /^{4>hosk(^dcji 


FiotrRB 7. — Phase relationships for harmonic pitching oscillation. 

G.(0 = J V^Jr)e- <" -> dr 

CAT)e-*“^dr+a„ain (9) 


= icoa{t') ^ 


Now, as in the previous example, let fi^(r) = (?/,^(d)— Pi(r) 
so that equation (9) becomes 

Gi^{t') = a{t')Gi^Jf,')-iwa{t') Pi(r)e-t-dr (10) 


Again, introduce the nondimensional parameters <p and k, 
separate equation (10) into its real and imaginary parts, 
and let <p approach infinity for subsonic speeds and <Xa for 
supersonic speeds. There results. 


R.P.(^^=Gi,J}i)—k ^ Fi{ip) sinkipd<p 

X= ro , Mo<l 
\=Ca, M„>1 

I.P.(^^=—k Fi{<p) cos k(pd<p 


( 11 ) 


> This quantity is, of course, the same lift force in phase with the pitching velocity which 
would occur alone had the wing been executing a steady turn (g constant). One of the chief 
advantages of the indicial response method, at least for supersonic speeds, is the case with 
which the relative importance of the various terms contributing to the total lift and moment 
can be assessed and the sources of the important contributions identified. 
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Notice in equations (11) for supersonic speeds that when the 
trigonometric terms are expanded for the low frequency 
case, as was done in the previous example, there appears an 
in-phase term of zero order in k, and an out-of-phase 

term of first order in k, —ika^ ° Fi{<p)dip. These, then, are 

the principal contributions to the lift forces for the low- 
frequency angle-of-attack variation. The phase relation- 
ships for this motion are shown graphically in figure 8. For 
this case, it is evident that the total-lift force can lag behind 
the angle of attack because of the negative out-of-phase 

contribution, —ika. f ° Fx{(p) cos k<pd<p. 


ika 



Figure 9. — Phase relationships for harmonic rotary oscillation. 


ika 



Figure 8. — Phase relationships for harmonic plunging oscillation. 

APPLICATION OF INDICIAL FUNCTIONS TO HARMONIC ROTARY 
OSCILLATIONS 

Finally, consider the case of harmonic rotary oscillations. 
Here, as previously mentioned, the normal velocity over the 
wing surface is composed of contributions from both the 
angular velocity and the instantaneous angle of attack, so 
that the complete expression in equation (3) must be em- 
ployed to obtain the total lift. However, for single-degree- 
of-freedom rotary oscillations, a equals 6 and a equals q, so 
that in this case the separate expressions given for the 
harmonically pitching wing (eqs. (7) and (8)) and the 
harmonically plunging wing (eq. (11)) can be combined to 
give the total lift for a wing executing harmonic rotary 
oscillations. Then, adding the results of equations (7), (8) 
and (11), the in-phase and out-of-phase components of the 
total lift become 

R.P. ^ kipdip-\-F ^ Fi{(p) cos ktpd<p 

(12a) 

X= to , 


It will be noted in figure 9 that the total-lift force can either 
lag behind or lead the angle of attack, depending on the 
relative magnitudes of the three terms comprising the 
out-of-phase lift. The total lift is shown lagging behind the 
angle of attack in figure 9, which situation, for axis positions 
ahead of the point of concentration of the total lift, gives rise 
to the possibility of the development of divergent rotary 
oscillations. 

Again, the complete frequency-dependent equations for 
the total lift of a wing in supersonic flow due to the rotary 
oscillation (eqs. (12), X=o-„) may be reduced to first order 
in k for the low-frequenc}^ case in the same manner as was 
described in the two previous examples to give 


i?,P.(^)=C'.„(.„) 

I.P. (^)=* [c,.S.<ra) - F, {<p) 


(13) 


For all three examples, the same procedure may, of course, 
be used to obtain the pitching-moment coefficient. Only 
the pitching-moment equations for the rotary oscillation are 
presented below, since the correspondence between the lift 
and moment equations is obvious. For the rotary oscilla- 
tion case, then. 


R.P.(££j=C^JF)-k £ FM sin kvdv+ 
k'^ ^ Fi{ip) cos kipd^p 

X=ra,M„<l 

X=(7a, Mo>l 

l.P. (X)~^: JV4(^) mvkipdv- 


K (14) 


/.P. ^=A:^C£jX)— ArJ* F 2 {<p)s\nkipd<p— ^ Fiifp) cos, k<p<F^ 

(12b) 

The phase relationships for the rotary oscillation may, of 
course, also be obtained by directly adding the results given 
in figures 7 and 8. The result of this addition is shown in 
figure 9. 




where, as previously. 


P3(^)=a„(x)-C'™„(^) 

and 



USE OF INDICIAL FUNCTION CONCEPT IN ANALYSIS OF UNSTEADY MOTIONS OF WINGS AND WING-TAIL COMBINATIONS 7 


Again, reducing equations (14) for supersonic speeds (X=o-a) 
to first order in frequency, there results 


IP- (%)=* 


The complete frequency-dependent equations for the lift 
and pitching-moment coefficients for the rotary-oscillation 
case (eqs. (12) and (14)) describe completely the aerody- 
namic forces and moments resulting from the single-degree-of- 
freedom pitching mode. For the purposes of the present 
discussion, however, it is sufficient to limit consideration to 
the simpler first order in frequency results of equations (13) 
and (15). The significance of the effect of the higher order 
terms on the out-of-phase pitching moment will be examined 
in a later section of this report. 

With regard to the first-order results for the supersonic 
lift and pitching-moment coefficients, it is instructive to 

note that the quantities j °Fi{ip)dip and [ in 

Jo Jo 

equations (13) and (15) are represented geometrically by the 
areas of the shaded portions of figure 10. 



Figure 10. — Geometric representation of the functionsj and 

r Ps(<t>)d<)>. 

0 


however, the anamolous result can be considered to be a 
consequence of assuming a two-dimensional flow, and there 
is reason to believe that the difficulty as the frequency 
approaches zero will not exist for finite-span wings. 

As has been mentioned previously, the use of the character- 
istic areas I Fi{tp)d(p SLJxd I will be shown to be of 

Jo Jo 

considerable value in estimating the damping-in-pitch 
characteristics of wings at supersonic speeds. For the two- 
dimensional wing at subsonic speeds, the singularity as 
k-^Q prevents the use of such a simplified approach without 
further study. However, rather than return to the use of 
the full frequency-dependent equations ((12) and (14)), the 
reduction of the equations for the out-of-phase lift and 
moment to first order in frequency will be made in such a 
manner as to preserve the significance of these areas. To 
accomplish this end, equation (12b) is reconsidered. It is 
evident that the first integral in equation (12b) may be 
discarded, since its contribution to the out-of-phase lift is of 
third order in frequency. The second integral is divided 
into two parts: 


J Fi{(p) cos k(pd<p= I 

o Jo 


Fx{<p) cos k(pd(p-\- 

J Fi{(fi) coa kipd<p (16) 


where tpi is chosen such that Fi(^i) is close to zero. The 
first integral in equation (16), being bounded, then causes no 
difficulty. Expanding the trigonometric term and retaining 


only the first term, in the expansion, there results 



which is the characteristic area out to the point <pi. Now 
for large values of ip, Fi((p) is approximated in reference 4 

by 




The manner in which these areas are affected by variations 
in Mach number, aspect ratio, plan-form shape, and fre- 
quency will be used as a guide in later sections of this report 
to determine the significance of these parameters. 

In the foregoing discussion, no mention has been made of 
reducing the complete equations for the lift and moment 
coefficients at subsonic speed (eqs. (12) and (14), X= qd) to 
first order in frequency as was done for the equations noted 
as applying at supersonic speed. It is evident that if the 
same procedure had been applied for subsonic speeds, the 

area corresponding to the term /; Fi(<p)dip can either be 

finite or can become infinitely large, depending on the 
manner in which the indicial lift function Ct„(¥>) approaches 
its steady-state asymptote as (p-^ ao . In the latter case, 
there exists the interesting anamoly of an infinite out-of-phase 
lift force as the frequency approaches zero. As can be seen 
from the results of reference 4, such will be the case for the 
two-dimensional wing. This result as the frequency ap- 
proaches zero is not peculiar to the indicial anaWsis alone, 
but has been pointed out by a number of authors using 
different approaches. As indicated by Miles in reference 10, 


where the values of a, and b are dependent on Mach 
number, and are given for M„=0, 0.5, and 0.8 in reference 4. 
Inserting this quantity in the second integral in equation 
(16), we have 



cos k<pdip=) 




cos kipdxp 
a-\-<p 



cos k<pd(p 


(17) 


Performing the indicated integrations in equation (17), 
there is obtained a term, from the second integral, 

and a term, — iiCi[k{a-'r (pi)] from the first integral. For 
small values of the argument, the cosine integral is approxi- 
mated by (see ref. 11), 

% 

Ci{k{a-\- vdi] ~ln[yk{a-{- <p-^] 


k{(i-\~ <p\) — >0 

where T is Euler’s constant, 1.78107. Then, through the 
first order in frequency, 

I.P. [(7^/ CO)- J V.(^)d^d-(?(^„ *)] (18a) 


8 

where 
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G{tpi, k)= 

The out-of-phase moment for pitching about the leading edge 
follows from the above development^ with the added result 
of reference 4 that for large values of <p, 


Fz{<p ) « 


EM 

4 


Then, through the first order, 

I.P. \_Grai.^)-^yz{<p)d<p-\ G{^uk)\ (18b) 

Thus, after fixing <pi, choosing a (small) value of k, and 
computing G(^i,k), the finite areas corresponding to the 

terms — I 'Fi (^)fZ^and — j ’i ^3 ((p) canbe assessed in the 
Jo Jo 

same manner as will be done for the supersonic case. The 
advantages of such a procedure will be evident later. 

CORRESPONDENCE BETWEEN INDICIAL LIFT AND MOMENT ANALYSIS 
AND OTHER METHODS 

Before proceeding further with applications of the indicial 
response ro.ethod, it is appropriate to discuss the relationship 
of this approach to other widely used methods. 

Following the fundamental papers of Bryan and Routh, 
which introduced the basic differential equations of motion 
of rigid bodies and their stability criteria, the historical 
development of the theory of longitudinal motions of an 
aircraft evolved separately in two fields of research: dynamic 
stability and flutter. Workers in the dynamic stability 
field soon found that the longitudinal oscillations of a rigid 
aircraft in flight were generall.y of small reduced frequency. 
On this basis, the constants due to the aerodynamic proper- 
ties of the airframe which appear in the differential equations 
of motion were considered to be independent of frequency. 
As a flrst approach to the problem of obtaining the necessarv 
aerodynamie coefficients analytically, the instantaneous 


normal velocity distribution at the surface of the airfoil was 
assumed to be constant with time. The aerodynamic forces 
and moments arising from the fixed boundary conditions 
were then calculated using steady-flow theory. Later, this 
assumption was realized to be an over-simplification for the 
case of wing-tail combinations and an additional term, essen- 
tially correct to the first order in frequency, was added which 
accounted for the lag in the tail pitching moment caused by 
the time required for the vorticity discharge from the wing 
to reach the tail (see ref. (12)) . Since at low speeds the pitch- 
ing moment of the tail far outweighs all other contributions, 
the results from steady-flow theory together with the term 
accounting for the vorticity lag satisfactorilj^ predicted the 
dynamic longitudinal motions of wing-tail combinations, 
and it was concluded that the major aerodynamic effects 
had been taken into account. In recent j’ears, however, 
numerous authors (in partieular. Miles; see, e. g., ref. 10) 
have pointed out that the above-mentioned theory overlooks 
important contributions to the aerodjmamic forces and 
moments which, though still within the first order in fre- 
quency approximation, arise from time-dependent boundary 
conditions and must be calculated from unsteady-flow theory. 
It has been shown by these authors that with proper inclusion 
in the equations of motion of these coefficients, the deteriora- 
tion of damping in the short-period mode actually occurring 
for aircraft flying at speeds near the speed of sound can be 
successfully predicted. The consequences of the assump- 
tions involved in the classical dynamic-stability theory will 
be more evident from a brief review of the equation of 
motion and boundary conditions for the single-degree-of- 
freedom rotary oscillations of a rigid wing fljdng at constant 
supersonic speed. At the very outset, the assumption is 
generally made that the aerodynamic reactions to the motion 
of the airframe depend onh^ on the angular position and 
angular velocity and not upon angular accelerations or higher 
time derivatives. The equation of motion for the change in 
pitching moment following a displacement from an equilib- 
rium position is then written in the form of a power series 
in a and a: 


I 

qoSc 








r -1 

(ac/2V„r , 

r n 

(d:C/2F„)- 1 


1 2! ^ 

LKt)J 

2! 


2! ' 


cross-derivatives and 
higher-order terms. 




(19) 


It should be remembered that for the rotary-oscillation case, 
the airfoil is subjected to changes in both angle of attack a 
and angular velocity q, and that these motions produce 
normal velocity patterns at the airfoil surface which are 
different in character. Thus, although for the single-degree- 
of-freedom case, a and q are equal, nevertheless their sep- 
arate effects must be considered and it is therefore necessary 


to include both 




and 




in equation (19). 


5(gc/2F„) b(acl2l\) 

Next, if it is assumed that the moments are linearR de- 
pendent on their respective variables, the second and higher- 
order terms in equation (19) may be discarded and the re- 


maining partial derivatives considered as constants for the 
given wing. There remains, therefore, a linear second-order 
system with constant coefficients. In order to calculate the 
coefficients (termed stability derivatives) theoreticalh^ it be- 
came necessary, for lack of more refined theoretical methods, 
to assume that the instantaneous normal velocity of the flow 
at the surface of the wing was fixed with respect to time. 


Thus, the partial derivative 


dC„ 


could be calculated 


d(qc/2V„) 

as the pitching moment due to a constant pitching rate, 
that is, Cm^(a„), while the derivative becomes the 

pitching moment due to a constant angle of attack, that is, 
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Cm„(<ra)- As a consequence of fixing the normal velocity 
pattern in time, however, it was necessary to assume that 


the derivative 


i>Cm 


was zero. There was therefore no 


b(acl2V„) 

possibility for this theory to predict the occurrence of dy- 
namic instability for a wing alone, since the only damping 
term remaining is Cm^{<ra), which is always stabihzing. When 
the resection of constant normal velocity with time is lifted, 
however, the assumption is then made that the stability 
derivatives in equation (19) may be calculated separately 
by fixing each of the mdependent variables a, a, and q in 
turn with respect to time. The derivatives Cm (O and 

bC 

Cm^ifa) thus remain unchanged, but the derivative ™ 


b{acl2V„) 


(or Cm^) can now be included and calculated as the pitching 
moment due to a constant vertical acceleration, a It 
should be emphasized that while Cm^((ra) and Cm^(<ra) may be 
calculated from steady-flow theory by virtue of the assumed 
invariance with time of the normal velocity pattern, Cm^ 
must be calculated from unsteady-flow theory since for con- 
stant a the angle of attack varies linearly with time, as does 
the normal velocity of the flow at the surface. 

It is clear that since the stability derivatives in equation 
(19) are assumed to be independent of the frequency, the 
result for the aerodynamic pitching-moment coefficient is 
thereby limited to one that is correct only to the first order 
in frequency. For the case of single-degree-of-freedom 
harmonic rotary oscillations, then, the in-phase and out-of- 
phase components of the total aerodynamic pitching-m.oro.ent 
coefficient, correct to the first order in frequency, become 


R.P. 

I.P. 




( 20 ) 


By comparison with the first-order in frequency result from 
the indicial response analysis (eq. (15)), it is evident that 

the two results are identical if the quantity ~ J" ” Pii‘p)d<p 

can be shown to be equivalent to To show this 

equivalence, consider a wing, initially in level steady flight, 
which is suddenly forced down with constant vertical accel- 
eration 6tV„. As seen in figure 11, the angle-of-attack vari- 



Figurb 11. — Flight path corresponding to sudden uniform normal 
acceleration a V„. 


ation in this case is a=at', where a is a constant. Then 
applying the counterpart of equation (3) for the pitching 
moment 

CmSr)<t'-r)dr 

Inserting a(t' — r)=a{t' — t], and performing the indicated 
differentiation, 

Cm{t’)=a£'CmSr)dr 

Now replace by Cm^{t')—Fi{T) and let t' be greater 

than f'. 

Then 

Cm{t') = a{t>) Cm^iO - a J^'“ F,ir) dr 

and nondimensionalizing, by replacing t' and f' by cV/2F„ 
and c<rj2V„, we have, 


C,n{‘p) = Oi{<p)C,nS^'^~'^ ^ " ^ii<P)dlp ‘P>(Ta ( 21 ) 

Thus, the pitching mom.ent proportional to the constant 
vertical acceleration parameter ac/2Vo, which is synono- 
m.ous with the definition of the stability derivative Cm^, is 
found to be equivalent to the pitching-m.om.ent contribu- 
tion due to a for the first order in frequency rotary-oscillation 
case.^ Therefore, the results of the indicial response method, 
when reduced to the first order in frequency for supersonic 
speeds, are identical to the results from the familiar first- 
order theory used in dynamic stability work.^ 

Workers in the field of flutter, who were concerned with 
frequencies m.any tim.es those encountered in dynamic- 
stability analyses, required theoretical information showing 
the behavior of the forces and moments as affected by the 
frequency of oscillation, and therefore discarded the first- 
order theory for more precise methods of analysis. One of 
the most useful of these has been the “oscillating potential” 
theory, which is based on solutions of the time-dependent 
linearized equation of com.pressible potential flow for the 
case of harmonic motions. The in-phase and out-of-phase 
lift and moments are thereby determined, generally as 
functions of powers of the reduced frequency, aspect ratio, 
Mach number, and position of the axis of rotation. The 
application of this method, which developed in this country 
primarily as a result of Theodorsen’s work for incompressible 
flow (ref. 13) has recently produced a number of useful 
papers covering a wide variety of wings at supersonic speeds 
(see, e. g., refs. 14 to 18). It has been shown by a number of 
authors, in particular, Garrick, in reference 7, that through 
the use of superposition methods the results for the aero- 
dynam.ic coefficients obtained from the oscillating potential 


- By the same procedure, the stability derivative Ct • (^a) can be shown to be equivalent 


to the term 




iip)d<p. 


3 Notice in equation (21) that if a® is replaced by infinity, the results apply to subsonic 


speeds. For the two-dimensional wing, the analogy between Cm- and 


-/■ 


Fi(<p)d<p then 


gives only the previously mentioned singularity at infinity as k-^0. If the area correspond- 

r co 

Fi(ip)d>p were finite, however, the analogy would be equally useful for subsonic 


as well as supersonic speeds. 
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theory are wholly compatible with those of the indicial 
response method applied to harmonic motions (eqs. (12) 
and (14)). 

Thus the indicial response method may be applied with 
equal vahdity to problems in both the fields of dynamic 
stability and flutter. In many cases, however, application 
of the specialized approaches mentioned above may result 
in greater economy of effort. The chief virtue of the indicial 
response method is that it represents a fundamental ap- 
proach to the problem of unsteady flow, and affords valuable 
insight into the physical nature of the aerodynamic phenomena 
taking place. 

PHYSICAL CONCEPTS RELATING TO THE INDICIAL LOADING 

It has been shown that for even small frequencies, the 
pitching moment of an airfoil in harmonic rotary motion 
can lag behind the angle of attack of the airfoil. The mag- 
nitude of the lag depends on the character of the indicial 
response to a step change in angle of attack. It is therefore 
of interest at this time to re-examine the physical nature of 
the flow that contributes this lag. 

Consider first the lift and moment at the instant the angle 
of attack changes, assuming that previous to time zero, the 
wing has been flying a level path at zero angle of attack. 
At t' = 0, the wing begins to sink, without pitching, with 
constant downward velocity aV„ while maintaining its 
forward velocity. The angle of attack therefore changes 
discontinuously from zero to a constant a. At the same 
instant, the step change in angle of attack causes the emission 
of a compression wave from each point on the lower surface 
of the wing and expansion waves from points on the upper 
surface. In the infinitesimal time during which the starting 
action occurs, each section of the wing experiences the same 
impulsive force, and by equating the im.pulse to the m.omen- 
tum transmitted to the mass of fluid affected by the starting 
waves, the starting lift coefficient can easily be derived as 
iajMo (see ref. 3). During the infinitesimal starting time, 
the pressure disturbances from, the edges of the airfoil, 
propagated at the speed of sound, travel an insignificant 
distance and do not influence the remainder of the airfoil. 
The lift coefficient is therefore independent of the wing plan 
form. This remarkably simple result for the starting lift 
coefficient, which is valid for both subsonic and supersonic 
speeds, is thus dependent solely on the flight Mach number. 
The starting pitching moment follows directly fi'om the 
above result, since by virtue of the uniformity of loading the 
aerodynamic center is located at the wing centroid of area. 

For values of time greater than zero, however, the situ- 
ation differs radically for the supersonic and subsonic speed 
ranges. Consider first the supersonic case. As time passes, 
the spherical sound waves emitted at t' = 0 grow in size 
with radius The wing, however, is moving forward at 
a faster rate than the rate of growth of the starting sound 
waves and thus begins to emerge from the influence of these 
waves. This is shown schematically in figure 12. (Note in 
fig. 12 that the wing moves away from a coordinate system 
that is fixed in space at the original position of the wing 
leading edge.) 


Leading edge at 



Figure 12. — Relation of wing position to starting sound waves for 
supersonic speed. 

At t' = 0, the starting waves just cover the wing and the 
loading is uniform as described previously. At = the 
starting waves have grown in radius and the wing has begun 
to emerge from their influence. On that portion of the wing 
which has emerged, region ® in figure 12, the loading has 
already reached its steady-state value. Notice that in this 
region the characteristic tip Mach cone has already formed. 
On the portion of the wing uninfluenced by the starting 
waves from the edges, region @ in figure 12, the loading is 
still uniform as at t' = 0. In regions ® and .© ,the loading is 
influenced by the starting waves from the leading and side 
edges, and in these regions is thus different from the loading 
in either region © or ®. As time increases still further, the 
uniform starting load quickly disappears as the sound waves 
from the leading edge grow in size and as the wing moves 
forward. Finally, at time the envelope of the starting 
waves from the leading edge is coincident with the trailing 
edge of the wing, and the steady-state loading corresponding 
to the new angle of attack a. has been completelj'' established 
over the wing. 

The above relationships can be shown more clearN for the 
entire time interval zero to <' for a two-dimensional wing 
by plotting as a function of time the position of the wing 
leading and trailing edges and position of the envelopes of 
the sound waves which emanate from the leading and trail- 
ing edges at <' = 0. Such a plot is shown in figure 13. 

It is clear that at t' = t\ the regions of the wing ©, ©, and 
© correspond to the same regions at t' = t[ for the wing 
shown in figure 12. For f' = 0 and in region © the loading 
is uniform and is given by 4a/M„. Note also in figure 13 
that disturbances created by the vorticity shed into the 
airfoil wake at <'==0 and all subsequent times cannot in- 
fluence the flow over the airfoil itself. For and in 

region © the wing has out-stripped the starting waves from 
the leading edge and has attained its steady-state loading. 
For t' <Ct'a the chordwise loading is composed of combinations 
of the loading in each of the three regions shown in figure 
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Figure 13. — Relation of wing position to starting sound waves for 
supersonic speed. 



Figure 14. — Indicial loading on two-dimensional wing at supersonic 
„ speed. 


14.^ Since the loading on the wing attains its final steady- 
state distribution at precisely the time when the wing has 
emerged entirely from the influence of the starting sound 
waves from the leading edge (or apex), the time to reach 
steady state may be easily calculated for any type of wing 
by means of the geometric relationships shown in figure 15. 
It may be easily verified that t', the time required for the 
wing to attain its steady-state loading, is given by either of 


< Tbe reader will note the similarity between figure 13 and sketches depleting the bound- 
ary conditions for three-dimensional wings in steady supersonic flow. Many researchers 
have pointed out the analogy and it has been used to calculate the pressure over a w ing 
i mpiilsively starting from rest (refs. 3 and 4). 



l 


Figure 15. — Geometric relationships between wing and starting sound 
wave at time t„. 


the following relationships; 


.1 Mo>('SC V-\ — J- 
CSC v) 0 


cos V, v< - 



( 22 ) 


The second of equations (22) applies to that range of Mach 
numbers for which the trace of the starting sound wave from 
the apex is not tangent to the trailing edge at Notice 

that for wings having straight or sweptforward trailing 
edges and straight or sweptback leading edges, equation (22) 
reduces to 


CoMo 

Vo{Mo-l) 


(23) 


In terms of the number of half M. A. C. lengths of travel, it 
thus appears that for wings having straight or sweptforward 
trailing edges and straight or sweptback leading edges, 


2Vo ,/ 2c„ Mo 

c (Mo-1) 


(24) 


whereas for wings having sweptback trailing edges and 
straight or sweptback leading edges. 


2c„ Mo , , - , 2co - ■X 

a-a=-^ rjTf > CSC v+-^ COS v,v<-^ 

c {Mo— CSC V) 6 2 


0-0= 2M„^ 




(25) 


^ 1 

< CSC V-\ — i- cos V 
~ b 


Now consider the subsonic case. Here the situation is 
more complicated in that, since the starting sound waves 
travel faster than the wing, the wing never escapes their 
influence. Furthermore, the vorticity shed by the wing at 
t'=0 can also influence the loading on the wing since the 
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disturbances created by the shed vorticity travel forward at 
a faster rate than the wing. For these reasons, the indicial 
loading at subsonic speeds approaches its steady-state dis- 
tribution asymptotically with time. The situation for sub- 
sonic speeds will be more clearly evident from examination 
of figure 16, which shows the relationship of the traces versus 
time of the leading and trailing edges of a two-dimensional 
wing flying at a subsonic Mach number to the traces of the 
envelopes of the starting sound waves. 



Figure 16. — Relation of wing position to starting sound waves for 
subsonic speed. 

Notice in figure 16 that the starting sound waves intersect 
the edges of the wing and that each intersection causes a new 
sound wave to be emitted, which in turn will intersect an 
edge. Furthermore, notice that the vorticity shed from the 
trailing edge at time zero can influence that portion of the 
wing behind the sound-wave trace labeled a-b. The influ- 
ence of each successive sound-wave reflection, however, is 
weaker than the last, and as the wing moves away from the 
starting vortices their influence diminishes, so that at <» 
the loading on the wing attains its steadj^-state distribution. 
The variation of the chordwise loading with time for the two- 
dimensional wing flying at a subsonic Mach number is shown 
in figure 17. Notice that for f>0, the chordwise loading is 
markedly different from the loading at supersonic speeds (fig. 
14). However, for t' — Q and in the region corresponding to 
region ® of figure 16, the loading is uniform and equal to 
AalMo, as in the supersonic case. 

DAMPING IN PITCH OF LOW-ASPECT-RATIO WINGS 

Previously (eq. (20)), it was shown that for single-degree- 
of-freedom, low-frequency, rotary oscillations of a wing, the 
principal parameter contributing to the damping of the 
motion is the damping coefficient This result, 

however, is not directly applicable to the case of an aircraft 
in flight, since generally additional damping is provided by 
virtue of the fact that the aircraft experiences harmonic ver- 
tical translatory oscillations as well as the rotary oscillations. 
It can be shown, however, that although the effect of the 
translatory oscillation is usually to increase the total damp- 
ing, nevertheless, the parameter of primary importance in 
determining the magnitude and duration of the motion re- 
mains the damping-in-pitch coefficient + The re- 
mainder of this section is therefore concerned with a study of 



the effect on this parameter of certain important variables. 
In particular, the effect of the position of the center of 
gravity, and the effects of Mach number, aspect ratio, plan- 
form shape, and frequency will be examined, principally by 
inspection of the indicial lift and moment responses to a 
change in angle of attack. 

EFFECT OF STATIC MARGIN 


From the previous discussion it will be remembered that 
at supersonic speeds the stability derivatives Ci,^(a-a) and 
^mi(va) were shown to be equivalent to the indicial lift and 


moment expressions. 


i' 


' Fi{>p)dip and 



Furthermore, it was shown that these were the contributions 
which could cause the total lift and moment during the short- 
period oscillation to lag behind the angle of attack. Hence, 
by inspection of equation (15), it is evident that since 
is always stabilizing, when is negative 

(corresponding to a statically stable condition) the possibility 
of d 3 mamic instability in the form of divergent rotary oscilla- 
tions arises when Cm^{(Ta), the shaded area in figure 10 (b), is 


larger than Cm„(o-J. Now since the normal velocity at the 
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surface of the wing due to the instantaneous angle of attack 
is constant over the wing, the lift derivative which arises 
from this boundary condition is independent of the axis 
position and Cm^ will therefore vary linearly with axis posi- 
tion. This variation is illustrated in figui-e 18, where AXo 
represents the distance of the center of gravity from the 
aerodynamic center, measured positive forward of the aero- 
dynamic center. 



Figure 18. — Variation of Cm and C„. with a.xis position. 

a a 


The parameter Cz,^, on the other hand, is a direct function of 
the axis position, since it arises from a normal velocity dis- 
tribution that varies directly as the distance from the axis. 
The moment coefficient Cm,^ will therefore vary as the second 
power of the axis position, and describes the parabolic shape 
shown also in figure 18. It is evident from figure 18 that 
the sum of Cm and Cm^ will be a minimum at some value of the 
static margin, and that the sign of at that point 

determines whether or not a region of axis positions will 
exist over which the wing can experience negatively damped 
rotary oscillations. These qualitative statements may be 
written explicitly by considering the equation for the damp- 
ing in pitch about an arbitrary axis, 

c.,+c.,= (C.,),+ (CO -X [(C..).+C.J-2 (xX C^. 

(26) 

where again Ar„ refers to the distance of the center of gravity 
from the aerodynamic center, and the subscripted terms are 
referred to an axis through the aerodynamic center. 

Taking the derivative of equatiou (26) with respect to 
Axofc, and setting the result equal to zero, there is obtained 
the axis location at which the damping in pitch is a minimum 




d(AXol'c) 






Y (27) 


'When equation (27) is inserted into equation (26), the mini- 
mum value of the damping in pitch is given as 


and hence, a region of instability wiU exist if 

[(C..).+ (C..)J+g^ [(Ci.),+aj'>0 (29) 

If equation (29) is greater than zero, the boundaries of the 
region of axis positions over which instability is possible are of 
course given by setting equation (26) equal to zero and 
solving for AXofc. 

1 4(7.„ ^VL 4(7.„ J+L 2(7.„ 

(30) 

Notice in both figure 18 and equation (30) that for a given 
Mach number there will be two axis positions at which the 
damping in pitch vanishes. Then if the above procedure 
is carried out for a series of Mach numbers, one may trace 
out a curve as shown in figure 19 which forms the locus of 



Figure 19. — Typical supersonic stability boundary curve. 


Mach numbers and axis positions at which the damping in 
pitch is zero. This locus thus delineates the regions of 
Mach number and axis position for which dynamic insta- 
bility is and is not theoretically possible. Such loci, cover- 
ing a wide variety of wings at supersonic speeds, have been 
presented by a number of authors. Watkins, for example, 
presents supersonic boundary curves for rectangular and 
triangular wings in references 15 and 16. At subsonic 
speeds. Miles’ reduction of Possio’s development to first 
order in frequency (ref. 10) can be used to form a stability 
boundary curve for the two-dimensional wing for a given 
(small) reduced frequency. Such a curve is presented for 
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the entire Mach number range in figure 20, where, here 
-\-Xo is the distauce of the axis of rotation behind the leading 
edge and k, the reduced frequency, is 0.011 for subsonic 
speeds and approaches zero for supersonic speeds. Notice 
in figure 20 that at both subsonic and supersonic speeds. 



Mach number, 


Figure 20. — Single-degree-of-freedom short-period pitching stability 
boundaries for the two-dimensional wing at subsonic and supersonic 
speeds. 


the range of Mach numbers over which dynamic instability 
is possible is largest for center-of-gravity positions forward 
of the aerodynamic center. Further, the largest range of 
axis positions over which dynamic instability is possible 
occurs near Mo=l. Both of these characteristics have 
been shown to be true as well for three-dimensional wings at 
supersonic speeds (see, e. g., refs. 15 to 18). 


EFFECT OF MACH NUMBER 


Next consider the effect of Mach number on the damping 
in pitch of a two-dimensional wing with axis at the leading 
edge. The variation with Mach number of the indicial 
pitching-moment response to a change in a will first be exam- 
ined, using the information given in the previous sections 
and the indicial curves given in reference 4. At supersonic 

speeds, the manner in which— J" F 3 (<p)d<p, the area corre- 
sponding to is affected can then be assessed and compared 
with Cm„ (da). At subsonic speeds, use is m.ade of equation 
(18b). It is evident in equation (18b) that by fixing k and 

choosing pi, such that the quantity ^ G{pi,k) is the same 


at each Mach number, one is free to compare finite areas 



on an equivalent basis. 


As has been mentioned previously, the starting lift, at 
any Mach number, is AaJMo and is concentrated at the 
midchord. At M„=0, therefore, there is an initial infinite 
pulse in the pitching moment about an axis coincident with 
the leading edge after which the indicial curve drops to 
7 t/ 4 and begins to grow asymptotically toward its steady- 
state value TV 12. At low subsonic Mach numbers other 
than zero, the initial pitching moment is less than infinite 
but very large, being 2/M„, and then falls before growing 
toward the steady value 7t/2/3. As the Mach number in- 
creases toward 1.0, the starting pitching moment falls while 


the asymptotic value grows, until at Mg=1.0 the indicial 
curve becomes unbounded in asymptotic moment. 

As seen in figure 21, the effect of increasing the Mach 
number at subsonic speeds is therefore to increase rapidly the 
area corresponding to the destabilizing moment contribu- 

X «>i 

Fs{p)dp. 



Figure 21. — Variation with subsonic Mach number of the indicial 
pitching-moment response (0) for the two-dimensional wing. 
Axis at leading edge. 


In figure 21, k was chosen to be 0.011 and the values of pi 

1 

were picked such that—- 6{pi,k) was -f-4.88 at each Mach 


number. In the following discussion, the damping moment 


— J* '^F3{ip)dip-\-4: 


.88 will be referred to as (7^^ for convenience. 


At supersonic speeds, the initial value of the pitching 
moment 2/M„ continues to drop with increasing Mach 
number, but here the steady-state pitching moment also 
begms to fall and at a faster rate than the starting moment, 
being 2//3. Even more important, as the Alach number 
increases, the number of half-chord lengths traveled to reach 
steady state decreases rapidly, being 22, for example, at 
Af„=l.l, as compared to 4 at M„=2. As seen in figure 22, 
the area representing 0^^, therefore shrinks rapidly with 
increasing supersonic Macb number and becomes relatively 
unimportant at Mach numbers greater than 2. The trend 
of with Mach number through the range 0<(Af„<!2 is 
more clearly evident in figure 23. It is seen that Cm^ is 
positive, or destabilizing, throughout the Mach number 
range and that its effect is most important at Mach numbers 
near 1.0. Also shown plotted for comparison in figure 23 
is the variation of Cm^ with Mach number. When the 
parameters are added, it is evident that the damping moment 
for the two-dimensional wing with axis at the 


leading edge is destabilizing in the Mach number range 
•0<M„<1-414. 


EFFECT OF ASPECT RATIO 


To illustrate the effect of aspect ratio, it is convenient to 
compare the supersonic damping-in-pitch characteristics 
of a group of triangular wings having subsonic leading edges. 
The wings are of equal area and differ only in aspect ratio. 
As was done previously, the indicial lift responses to a change 
in angle of attack will first be examined. The effect of aspect 
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00, Mo =1.0 



Figube 22. — Variation with supersonic Mach number of the indicial 
pitching-moment response for the two-dimensional wing. 

Axis at leading edge. 



Figube 23. — Variation with Mach number of the damping in pitch of 
the two-dimensional wing. Axis at leading edge. 


ratio on the characteristic area l epi’esenting can then be 
assessed. 

As has been mentioned previously, the starting lift 
coefficient after a step change in a is independent of aspect 
ratio and is therefore equal to for each wing. The 

parameter Ua, the number of half M. A. C. lengths required 
to reach steady state, is also the same for each wing, being 
a function only of Mach number. The steady-state lift 
coefficient, on the other hand, is a function of aspect ratio 



Figube 24. — Effect of aspect ratio on the indicial lift response 
of triangular wings having subsonic leading edges. 


and decreases as the aspect ratio is reduced. Thus, as shown 
schematically in figure 24, as the aspect ratio becomes 
smaller, the characteristic area representing 0^^ decreases 
rapidly.® For the wing of smallest aspect ratio, Cl^ may be 
positive since the area below the steady-state lift coefficient 
is more than compensated for by the area above. It is 
evident, therefore, that a reduction in aspect ratio has a 
highly stabilizing effect on the damping in pitch, since for 
positive values of the static margin, the development of a 
destabilizing damping moment is possible only when is 
negative. This result is shown in figure 25 where, for an 
axis of rotation located at 0.25 c and M„=\.2, the damping 
parameters are presented as functions of aspect ratio. Since 
for triangular wings the lift due to a. is coneentrated at 
5/8 c, Cmg, is equal to —3/8 Cl^. The variation of Cm^ with 
aspect ratio shown in figure 25 then follows directly from the 
trend of shown in figure 24. Also plotted in figure 25 
is the variation of Cm^ with aspect ratio (ref. 19). 



Figube 26. — Variation with aspect ratio of the damping-in-piteh 
coefficients of triangular wings. Axis at 0.25c and Mach number 

1 . 2 . 

It is apparent that although Cm^ becomes more stabilizing 
with increasing aspect ratio, the destabilizing effect of C'm^ 
predominates, and the trend of the net damping mom.ent is 
seen to become highly destabilizing as the aspect ratio is 
increased. 

By the same reasoning, the variation with aspect ratio 
of the damping moments of other types of wings can be 
shown to be similar (see, e. g., refs. 10, 20, 21, and 22). A 
notable exception, however, is the triangular wing with 
supersonic leading edges, whose damping in pitch has been 
shown to be independent of aspect ratio (see refs. 17 and 18). 
This characteristic may be anticipated from a study of the 
indicial response curves, since not only are the initial pitching 
moment Cm^ (0) and the half M. A. C. lengths traveled to 
reach steady state (uo) independent of aspect ratio, but, 
unlilre the subsonic-edged triangular wing, the steady-state 
pitching moment Cm^ is also independent of aspect ratio, 

2 Theoretical indicial curves have not yet been calculated for the triangular wing with sub- 
sonic leading edges. The curves drawn in figure 24 are estimates of the true shapes, and are 
intended only to indicate the trend of the characteristic area with aspect ratio. The exact 
variation of Cl • aspect ratio can be computed from the results of reference 19. 
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4 / "\ 

being — - Inspection of the results of reference 4 

then reveals that the indicial variation Cm^ (<p) between zero 
and (Ta and the steady-state parameter Cm^ (va) are likewise 
independent of aspect ratio. 

EFFECT OF PLAN-FOKM SHAPE 

Next, consider the effect of plan-form shape on the super- 
sonic damping in pitch of a group of wings having the same 
aspect ratio. For this comparison, three wings of aspect 
ratio 3 are chosen, having the triangular, swept, and rec- 
tangular plan forms shown in figure 26. 



Figure 26. — Wings used in study of effect of plan-form shape. 


As has been shown in the section entitled “Effect of Static 
Margin,” the tendencies of the wings toward dynamic 
instability can be compared comprehensively by plotting 
their stability boundaries. For this comparison, then, use 
is made of equation (30). The stability derivatives which 
appear in equation (30) were computed from the theoretical 
results of references 10, 19, 21, and 22. Results of these 
calculations are shown in figure 27 wherein the stability 
boundaries for the three wings are shown as a function of 
axis position and Mach number. (Note that the axis posi- 
tion for each wing is measured as the distance from the 
leading edge of the M. A. C. of the wing, and that the 
dimensions are nondimensionalized on an equivalent basis 
by referring them to the M. A. C. of the triangular wing.) 
It is clear from inspection of figure 27 that at any Mach 
number the triangular wing has the smallest range of axis 
positions over which dynamic instability is possible and the 
rectangular wing, the largest. 

The differences in the damping characteristics of the 
triangular and rectangular wings will be more clearly under- 
stood by a qualitative study of their indicial lift response 
for a Mach number of 1.2, and an examination of the dis- 
tribution of loading due to a for the two wings. Consider 
first the indicial lift responses. 

Again, the starting lift coefScient is independent of plan- 
form shape and is 4a/Mo for each wing. For the rectangular 
wing, the lift drops abruptly after time zero due to the loss 
in lift in the regions of the wing influenced by the formation 
of the tip Mach cones and the starting waves from the side 
edges (see fig. 12). Then, as the wing begins to emerge from 
the influence of the starting waves, the lift begins to recover, 
rises toward its steady-state value (given by eq. 6.3-2 of ref. 
10), and attains this value after 12 half -chord lengths of 
travel (eq. (24)). The variation is shown in reference 4 and 
is redrawn in figure 28. As mentioned previously, theoretical 
indicial lift results have not yet been developed for triangular 
wings having subsonic leading edges. However, the varia- 
tion shown in figure 28 is considered to be a reasonable esti- 



Moch num ber. Mg 


Figure 27. — Comparison of single-degree-of-freedom short-period 
pitching stability boundaries at supersonic speeds for three wings of 
aspect ratio 3, having triangular, swept, and rectangular planforms. 


/}=3 





Figure 28. — Comparison of indicial lift responses at Mach number 1.2 
of rectangular and triangular wings of aspect ratio 3. 

mate of the true shape, being based on knowledge of the 
steady-state lift (ref. 19), the time to reach steady state (eq. 
(24)), and the assumption that the shape of the variation 
would be similar to that of the wide triangular wing (ref. 4) . 
The curve was adjusted within the known end points until 
the area corresponding to agreed with that given for 
this parameter in reference 19. It is evident from examina- 
tion of figure 28 that because of the initial loss in lift and the 
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larger steady-state lift for the rectangular wing, for 
this wing is significantly more negative than that for the 
triangular wing. Next, it is shown in references 10 and 20 
that with the exception of regions influenced by tip Mach 
cones, the loading due to a (for ip o-„) for wings having 
sweptback leading edges increases linearly from zero along 
rays from the apex; whereas for rectangular wings the load- 
ing due to a -increases linearly from zero along chord lines. 
These characteristics place the center of loading due to ct 
at % Co for the triangular wing and approximately % c<, for 
the rectangular wing.® Then for an axis of rotation passing 
through the aerodynamic centers of the wings, the moment 
arm for the lift due to a for the triangular wing is K 2 or 
% c ; whereas for the rectangular wing it is approximately % 
CoJ Thus, not only is the negative out-of-phase lift contri- 
bution (7x,^ for the rectangular wing significantly larger than 
that for the triangular wing, but the destabilizing damping 

moment — larger yet, due to the larger moment arm. 

Calculations for the steady pitching parameter for an 
axis through the aerodynamic center (refs. 10 and 19) then 
reveal that (7m„ for the triangular wing is more negative than 
for the rectangular wing. The net result is therefore a con- 
siderably larger damping moment for the triangular wing 
than for the rectangular wing. The result of this compari- 
son, however, should not be interpreted as a recommendation 
that the triangular rather than the rectangular wing be used 
on aircraft from a dynamic stability standpoint. To obtain 
adequate static stability, the rectangular wing would gen- 
erally be employed in combination with a tail surface; whereas 
the triangular plan form m.ay be sufficiently airworthy with- 
out the use of a tail. The addition of a tail surface in effect 
reduces the aspect ratio of the rectangular wing, which re- 
duction, as noted previously, has a highly stabilizing effect 
on the damping in pitch. The tailless triangular wing may 
therefore experience more difficulty at Mach numbers near 
1.0 than a rectangular wing-tail combination. 


EFFECT OF FREQUENCY 


The previous discussion has been restricted to the analysis 
of a harmonic motion that is of vanishingly small frequency. 
This limitation arose as a consequence of discarding all but 
first order in frequency terms in the expansions of equations 
(12) and (14). The question arises: When the frequency 
can no longer be considered small, what effect has the 
frequency on the damping in pitch? 

Previously, the trigonometric terms in equation (14) were 
expanded and, assuming k to be very small, terms of 
order k^ and higher were eliminated. It was then found 
that the loss in damping from that provided by the steady 
damping parameter was associated with the de- 

stabilizing moment contribution corresponding to the term 


-j; 


Fi{<p)d<p. 


Now, however, we discard the restriction 


of small k and perform graphically the integrations evident 


* Due to the influence of the tip Mach cones, the center of loading due to a is shifted for- 
ward somewhat from the position It has for the two-dimensional wing. Calculations for the 
A=Z rectangular wing at Afo=1.2 show that the center of loading Is at 0.605c«. 

» Again, due to the influence of tip Mach cones, the aerodynamic center is shifted forward 
from c* to 0.443 c«. The moment arm is therefore 0.162 c«. 



Figure 29. — Illustration of graphical integration procedure for finding 
effect of frequency at supersonic speed. 



Figure 30. — Effect of frequency on the single-degree-of-freedom short- 
period pitching stability boundaries at supersonic speeds for a 
triangular wing of aspect ratio 4. 


in equation (14) for supersonic speeds for several values of k. 
The procedure is indicated in figure 29. It is apparent 
from figure 29 that the effect of increasing k is to 
reduce the area corresponding to the destabdizing moment 

X Oa 

F 3 (<p)cosk<pd<p. There appears another 
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destabilizing contribution, Fi{<p)smk<pd(p, but quite 

evidently its effect is small compared to tbe reduction in 
F^{<p)cofik<pd<p. Notice further in figure 29 


the term 


-j; 


that the effect of increasing k becomes of m.arked impor- 
tance when the half-period of oscillation is the sam.e order of 
magnitude as the time for the indicial response to reach 
steady state. As shown in figure 29 for the frequency k^, 
the area h then begins to subtract from a, so that the de- 


stabilizing contribution 




Fi{tp)cosk<pd((> can be very 



Figure 31. — Effect of frequency on the single-degree-'of-freedom short- 
period pitching stability boundaries at supersonic speeds for a 
rectangular wing of aspect ratio 4. 


small. We may therefore expect that increasing the fre- 
quency of oscillation has a stabilizing effect on the damping 
in pitch. This conjecture is substantiated in figures 30 
and 31, where the supersonic stability boundary curves 
for aspect-ratio-4 triangular and rectangular wings are 
shown plotted for various reduced frequencies. These 
curves were obtained from calculations based on the results 
of references 15 and 16. Notice that for both wings the 
region of possible instability is diminished as the frequency 
is increased. 

From the results of the analysis for supersonic speeds, we 
may further expect that the stabilizing effect of increasing 
the frequency Avill be of even m.ore importance at subsonic 
speeds, for here the indicial variation F<i{(p) dies out at 




Figure 32. — Illustration of graphical integration procedure for finding 
effect of frequency at subsonic speed. 



Figure 33. — Variation with reduced frequency of the single-degree-of- 
freedom rotary damping-moment coefficient for the two-dimensional 
wing at Mach numbers 0, 0.50, and 0.80. Axis at leading edge. 


infinity. The half-period of oscillation is therefore always 
smaller than the time to reach steady state. The situation 

for the term — J Fi((p)cos ktpdip is shown in figure 32. 
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It is evident in figure 32 that the destabilizing moment 

CO 

Fz{ip)cosk<pd(p diminishes rapidly as the frequency is 

increased. The effect of this reduction on the damping in 
pitch can be illustrated by plotting the subsonic damping- 
moment coefficient against reduced frequency for the two- 
dimensional wing (with axis at the leading edge) for Mach 
numbers 0, 0.5, and 0.8. The results, which were ob- 
tained from reference 4, are shown in figure 33. It is seen 
in figure 33 that the large destabilizing effects of the mo- 


ment contributions — I F^{tp)cos,k<pd(p and —k I i^ 4 (^)sin 

Jo Jo 

ktpdip are confined to a relatively narrow range of reduced 
frequencies. Notice further in figure 33 that the range of 
frequencies for which instability is possible is small at 
M„—Q (0<i<0.04) and grows with increasing Mach 
number. This is believed to be the primary reason why 
unsteady lift effects were found to be unimportant at low 
speeds but are of great importance at speeds near the 
sonic speed. 


PART II— WING -TAIL COMBINATIONS 


In the historical development of the field of dynamic 
stability, there has been general agreement among researchers 
that for conventional tail-aft aircraft, the primary source of 
pitch damping is that provided by the tail. In fact, several 
early investigators concluded that the tail contribution was 
the only one of consequence and, hence, that the contri- 
butions of the wing and fuselage to the total damping in 
pitch could be neglected with little error (c. f., ref. (12)). 
It was also recognized that a major portion of the lift de- 
veloped by the tail arose from the influence on the tail of the 
vorticity shed bj^ the forward wing, and, in particular, that 
a finite time was required for this vorticit 3 ^ to reach the tail, 
thereby creating a lag in the development of tail lift behind 
the motion of the aircraft. These considerations led to the 
formulation of a simple tail damping factor which took into 
account the time lag in vorticity in terms of the configuration 
geometry and the effective angle of attack at the tail, and 
which alone proved to be adequate for predicting the 
dynamic-stability characteristics of low-speed aircraft. 

It has already been shown in previous sections of this re- 
port that for higher speeds, the damping moment provided 
by the wing surface itself (and, analogously, the tail surface) 
can be highly significant. In subsequent sections the- im- 
portance of the influence of the vorticity shed by the forward 
wing on the damping moment provided by the tail is re- 
examined. Again, the indicial-function concept is used, 
and proves to be of valuable assistance in illustrating the 
physical nature of the problem. 

GENERAL CONSIDERATIONS 

THE FOUK COMPONENTS 

For each of the three harmonic motions considered earlier, 
the full frequency-dependent equations for the lift and 
moment developed for a single wing and their first-order in 
frequency counterparts are still applicable for a wing-tail 
combination, with the following reinterpretation of the in- 
dicial functions C'z.Jv’), C'»„(?>), and (7„j(^): each 

of the above functions must now be considered to rep- 
resent the combined responses of the wing and tail plus 
interference effects between the wing and tail. Obviously, 
however, since the theory is linear, the combined function is 
equal to the sum of its components. In the general case, 
four such components for each of the indicial functions may 
be enumerated. These are listed below for the indicial 
function {ip), that is, for the response in lift of the wing- 
tad combination to a step change in angle of attack of the 
combination. 


(1) The response in lift of the forward wing to a step 
change in angle of attack of the forward wing, the tail being 
at zero angle of attack. 

(2) The response in lift of the tail to a step change in 
angle of attack of the tail, the forward wing being at zero 
angle of attack. 

(3) The response in lift of the tail to a step change in angle 
of attack of the forward wing, thq tail being at zero angle of 
attack. 

(4) The response in lift of the forward wing to a step 
change in angle of attack of the tail, the forward wing being 
at zero angle of attack. 

The situation is entirely analogous for the other three in- 
dicial functions. It should be noted, however, that the 
components of the indicial functions due to pitching velocit 3 ^, 
Gh^iv) and G^J^tp), must all be referred to the same axis 
of rotation. In the following discussion regarding the four 
components, the use of the indicial lift response to a step 
change in angle of attack for illustrative purposes is continued. 
It will be understood that the remarks apph^ as well to the 
other indicial functions. 

For supersonic speeds, the component (1) is exactly 
equivalent to the response in lift of an isolated wing, and 
may be calculated on that basis, utilizing the theoretical 
results noted previously as applying in the supersonic speed 
range (refs. 3, 4, 5, and 6). Likewise, component (2) may 
be computed on the same basis. Furthermore, since at 
supersonic speeds, the disturbances created by the tail can 
never overtake the wing, component (4) is zero. There 
remains to consider, therefore, only component (3), for 
which-there are as yet no theoretical solutions. 

For subsonic speeds, as usual, the situation is more com- 
plicated. The component (1) is no longer rigorously equiv- 
alent to the lift response of an isolated wing, since the 
disturbances created by the wing by its own motion cause 
disturbances to be created at the tail which travel forward 
at the speed of sound, overtake the wing and, in turn, in- 
fluence the development of lift on the wing. Likewise, the 
reverse situation exists for component (2), so that the re- 
sponse in lift of the tad. is no longer that of an isolated tad. 
Further, component (4) is not zero, again by virtue of the 
fact that disturbances propagated at the tad are able to 
overtake the wing. lor subsonic speeds, therefore, all. 
four components require new solutions. It may be expected, 
however, that the secondary coupling effects between the 
wing and tad, as exemplified by component (4), are of small 
magnitude, and may usually be neglected. If it is assumed 
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that component (4) is zero, it is within the same order of 
approximation to assume that the coupling effects on com- 
ponents (1) and (2) are also zero, so that the subsonic 
nterference problem is then reduced to the same form as that 
of the supersonic problem. In the interest of brevity, 
however, in what follows the subsonic case will not be 
considered, except to note that the general procedures to 
be undertaken for obtaining numerical solutions for the two- 
dimensional supersonic case would be the same for the 
subsonic case, albeit, considerably more complex.® 

In the next section, therefore, an heuristic discussion is 
devoted to the remaining new function, component (3) for 
supersonic speeds; that is, the response in lift of the tail to 
a step change in angle of attack of the wing. 

BOUNDARY CONDITIONS AT THE TAIL CORRESPONDING TO A STEP 
CHANGE IN ANGLE OF ATTACK OF THE WING 

The boundary conditions corresponding to component (3) 
may be illustrated clearly by plotting as functions of time 
the traces of the wing and tail leading and trailing edges, 
and the traces of the sound waves emitted at the wing 
leading and trailing edge at the start of the motion. No 
sound waves are propagated at the tail at i' = 0 since, under 
the conditions of the problem, the tail remains at zero angle 
of attack throughout the motion. The plot is shown in 
figure 34.® 



Figure 34. — Boundary conditions at tail corresponding to step change 
in angle of attack of wing. 


It is clear from examination of figure 34 that at the start of 
the motion and until time t[, the disturbances propagated 
at the wing as a consequence of the change in angle of attack 
have not signaled their presence to the tail. During this 
time interval, obviously no lift is developed by the tail. For 
time greater than t[, the tail begins to penetrate the upwash 
field created by the vorticity shed by the forward wing, and 
thereby begins to develop positive lift on that portion of its 
surface which has penetrated the field. As time increases 

« It is to be noted that Jones and Fehlner have treated the incompressible case (ref. 23) 
by a method very similar to that described in a subsequent section of this report entitled, 
“Application of Methods of Oust Analysis.” 

9 It is implied In figure 34 that the wing and tail are coplanar or nearly coplanar surfaces. 
This restriction, while not a necessary one, is advisable from the standpoint of simplicity, 
and will hereinafter presumed to be in force. 


beyond the tail penetrates the field of shed vorticity 
itself and experiences combinations of both upwash and 
downwash on its surface. As the tail penetrates further into 
the field of vorticity, the region of downwash will predom- 
inate, so that it may be expected that the lift developed by 
the tail will eventually change from positive to negative. At 
time t's the wing has completed its build-up in lift which 
thereafter remains constant with time. The build-up in lift 
at the tail, however, continues until f', beyond which time 
the pressures on the tail are no longer influenced by the 
disturbances created by the wing at the start of the motion. 
For time greater than i', therefore, the lift on the tail 
remains constant with time and is equal in magnitude to the 
steady-state interference lift on a tail immersed in the non- 
uniform but steady downwash field created by the 
forward wing. 

From the above qualitative considerations, one may an- 
ticipate that the variation with time of the interference lift 
of the tail in response to a step change in angle of attack of 
the wing will resemble the variation shown in figure 35. 



Figure 35. — Qualitative prediction of variation with time of tail lift 
response to step change in angle of attack of wing. 


For wings and tails having straight trailing edges, the values 
of time t\ and f' in figure 35 may be computed exactly from 
the geometric relationships shown in figure 34. They are, 

Mo^ 

Vo{\+Mo) 

Mo (ciF+Cr+f) 

F„(M„-l) ’ 

Furthermore, theoretical solutions exist for the steady-state 
downwash behind wings of a wide variety of plan forms 
from which the final value of the interference lift (Q 
may be computed. It remains, therefore, to fix with a 
greater degree of certainty the variation (<') between 
t'l and f'. In view of the paucity of information concerning 
the nature of this variation, the succeeding six sections 
of this paper are devoted to the development of theoretical 
techniques from which results for certain representative 
cases can be obtained. 

THE TWO-DIMENSIONAL CASE 

BOUNDARY CONDITIONS 

Consider first the problem of calculating the growth in 
lift at supersonic speed on a tail surface of infinite span in 
response to a step change in angle of attack of a coplanar 
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forward surface, which is also of infinite span. This problem 
in two dimensions is chosen as a starting point since, here, 
use can be made of a valuable theorem; namely, that any 
unsteady-flow problem involving only two spatial dimensions 
can be related to and solved as an analogous problem in a 
supersonic three-dimensional steady flow. The basis for 
this statement can be made clear by a comparison of the 
governing linearized perturbation velocity potential equa- 
tions refefied to fixed-space axes for unsteady two-dimen- 
sional flow and steady three-dimensional supersonic flow. 
The former is written as, 

( 32 ) 

Q'o 

and is transformed to the canonical form of the wave equa- 
tion by the substitution t=aj,' : 

CSS) 

For the steady three-dimensional supersonic case, the 
potential equation is 

(34) 

It is then apparent that simply by replacing x, y, and jS in 
equation (34) by t, x, and 1, respectively, equation (34) can 
be made identical to equation (33). Thus, a solution to 
equation (34) for a Mach number of ^ (j3=l) can likewise 
be taken as a solution to equation (33) with proper reinter- 
pretation of the variables x, y, z. A complete development 
of this analogy can be found in references 3 and 4 along with 
several applications to problems of isolated wings. It 
suffices here simply to state that the theorem may be applied 
with equal validity to problems of wings in combination. 
In the light of this result, figure 34, which depicts the bound- 
ary conditions for the two-dimensional unsteady flow 
problem, takes on the following interpretation for the 
analogous three-dimensional steady-flow problem; the time 
axis becomes a spatial coordinate alined with the free-stream 
direction; the x axis becomes the spanwise coordinate, normal 
to the free-stream direction. The z axis is normal to the 
xt' plane and again is taken as positive downward. The 
traces of the wing and tail leading and trading edges and the 
projections of their chords on the x axis become in the steady- 
flow problem the projections onto the xt' plane of the physi- 
cal outlines of a pair of semi-infinite sweptback wings with 
tips normal to the free-stream direction. In the unsteady- 
flow problem the forward wing is caused to attain at t'=0 
an angle of attack a which remains constant for all time 
thereafter. This boundary condition is interpreted in the 
analogous steady-flow problem to mean that a uniform 
normal velocity of magnitude V„a exists everywhere over 
the surface of the forward wing. The traces of the sound 
waves which are emitted from this wing at f' = 0 as a conse- 
quence of the change in angle of attack become in the steady 
supersonic problem the characteristic Mach lines which 
emanate from the foremost extremities of the wing. The 
rear wing, being at zero angle of attack in the unsteady 
problem, remains so in the analogous steady problem, and 
thus no Mach lines appear at the extremities of this wing in 


the steady-flow analogy. The procedure for calculating the 
lift at the tail in response to the change in angle of attack of 
the wing is then theoretically clear. By any of the known 
methods available in steady three-dimensional supersonic 
flow theory one computes in the region occupied by the tail 
the perturbation-velocity field caused by the presence of the 
forward lifting wing. In the plane of the tail surface, this 
field contributes a normal-velocity fiow component which is 
canceled by the development of loading at the taU. A 
“spanwise” integration of the loading over the tail then 
yields the desired result — ■ tail lift as a function of time. 

SIMPLIFIED PBOBLf;M 

The procedure outlined above in principle constitutes a 
straightforward exercise in the application of the three- 
dimensional linearized supersonic steady-fiow theory. Un- 
fortunately, in practice, the multiple integrations that arise 
in the course of the computations are exceedingly complex, 
and have thus far defied solution. These computations were 
therefore discontinued in favor of a simpler approach which, 
in itself, retains the essential nature of the problem and has 
the additional advantage of lending itself readily to extension 
and further refinement. The simplification consists of 
replacing the forward wing by a single two-dimensional 
vortex which, just as the wing, starts impulsively from rest 
at f'=0 and moves forward in the 2=0 plane with constant 
supersonic velocity V„. In addition, in satisfaction of the 
condition that a vortex cannot end in a fluid, there must be 
placed at the point of departure of the moving vortex, a 
stationary or starting vortex of equal and opposite strength.''* 
It is well-known that the boundary conditions for a lifting 
surface may be satisfied by an appropriate distribution of 
vortices; thus the result to be obtained from this simplifica- 
tion m.ay be viewed as a fundamental one in the sense that 
the solution to the original problem may in principle be 
recovered by an appropriate superposition of elementary 
solutions. A procedure for accomplishing this task is de- 
veloped in a later section of this paper. 

The boundary conditions for the simplified problem are 
illustrated in figure 36. Here, “true time” t' has been re- 
placed by its spatial analog t, the relationship being t=aol' . 
The boundary condition for a closed vortex system is that the 



Figure 36. — Boundary conditions for simplified problem. 


10 The extremities of the movlog and starting vortices must, in general, be connected by 
two trailing vortices to form the usual closed loop. However, for the two-dimensional prob- 
lem, it is supposed that the trailing vortices are so far away from the centers of the moving 
and starting vortices that their effect on the flow there may be neglected. 
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jump in potential is a constant within the loop formed by the 
moving, starting, and trailing vortices. For the two- 
dimensional case this condition translates in figure 36 to the 
stipulation that the potential difference A$ across the 2=0 
plane be a constant within the shaded region bounded by the 
traces of the moving and starting vortices. The problem 
of calculating the induced velocity flow field associated with 
the vortex-pair may now be related to the analogous three- 
dimensional steady-flow problem; namely, that of calculating 
the induced flow field due to the presence in a supersonic 
stream of a yawed triangular wing, one edge of which is 
parallel to the free stream, the other being swept ahead of the 
Mach cone from the apex, and whose surface is warped in 
such a way that the potential difference across it is a constant. 
This problem has been solved by Lomax, Heaslet, and Fuller 
in reference 24, wherein it is found that the induced vertical 
velocity in the 2=0 plane for supersonic speeds is given by 
the relation. 

Thus, the problem of calculating the lift at the tail in response 
to the step change in angle of attack of the forward wing has 
now been reduced to one which can be stated in terms of the 
steady-state analogy as follows: Find the section lift as a 
function of t on the semi-infinite sweptback plan form shown 
in figure 36 which rests at zero angle of attack in a supersonic 
stream and which is subjected to a vertical-flow distribution 
given by equation (35) on that portion of its surface lying 
within the region —t<x<t. 

METHOD OF SOLUTION 

The problem stated above may be solved by any of a 
number of methods available in steady supersonic flow 
theory. One that immediately suggests itself, however, by 
the conical nature of equation (35) is the method of super- 
position of elem:entary sectors. If in equation (35) the 
substitution is made, a—xit, we get, 

(36) 

The quantity a of course, defines a ray emanating from the 
origin of the xt coordinate system, and by equation (36) it is 
seen that along a given ray, the normal velocity varies 
inversely as x. Then if the normal-velocity distribution is 
decomposed into infinitesimal steps of magnitude dw along 
rays from the origin, it can be specified that within each semi- 
infinite sector of the tail formed by the tail leading and trail- 
ing edges and a ray from the origin, the normal velocity 



sector, given by, 

It is then a relatively simple problem to solve for the loading 
on a single sector caused by the normal- velocity distribution 


and thence by the superposition of elementary solutions, to 
find the loading and integrated loading caused by the given 
normal-velocity distribution. This superposition procedure 
is briefly outlined in the following paragraph. 

The boundary conditions for an elementary sector formed 
by an arbitrary ray a are illustrated in figure 37 . 



Figure 37. — Boundary conditions for elementary sector. 


Within the shaded region shown in the figure the normal 
velocity varies as elsewhere on the plan form it is zero. 

At? 

Presuming now that the loading coeflficient — {x,t,a) at the 

2o 

point x,t for the elementary sector has been calculated, it is 
desired to sum the effects of all sectors that influence the load- 
ing at the fixed point. A contribution to the loading at the 
point x,t will be made by each sector within whose zone of in- 
fluence the point lies. As seen on figure 37 this region ex- 
tends from the sector defined by the ray a=l to that sector 
whose trailing Mach line from the apex of the sector passes 
through the point. On figure 37 this sector is defined by the 
ray a=d{x,t)- The total loading at x,t is then expressed by 
the formula 

^ (r,f)= r ^ (r,<,l) (38) 

g_o j 1 g.0 gto 

At? 

where \j/ (a) is given by equation (37) and — (x,t,a) is the 

loading coefficient at a point on the elementary sector due 

to a normal-velocity distribution — The second term in 

equation (38) accounts for the effect of a possible jump in w 
at the ray a=l. Howevei-, in the present problem, tf'(l) 
is zero so that this term disappears. The section lift at the 
station t is then found simply by an integration of the loading 
in the x direction, 

Ci(t)=- f — (x,t)dx=- f f ^'(a) — {x,t,a)dadx (39) 
There remains to discuss the calculation for the sector 
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loading — (x,t,a) at the fixed point x,t. This calculation, 

however, involving as it does a wing having all supersonic 
edges, is straightforward and does not require elaboration. 
An interesting sidelight arising as a result of this work was 
the discovery that the loading is zero at a point on the sector 
whose forecone includes only the sector leading edge and 
intersects the line a"=0. 

’ The section lift was found to vary in form in each of the 
five time intervals designated in figure 38. 



Figure 38. — Regions used in analysis of lift response of two-dimen- 
sional tail. 


Unfortunately, difficulties of integration prevent the presen- 
tation of solutions in closed form for the section lift in all hut 
the first time interval. It was always possible, however, to 
reduce the integrals to single quadratures, so that numerical 
or graphical methods can be employed to obtain results for 
specific cases. The results in integral form for the section 
lift in the intervals 1 to 5 are somewhat lengthy and therefore 
are presented in Appendix A. 

DISCUSSION 

In order to illustrate the nature of the result, the variation 
of section lift at the tail with “time” t has been computed 
for the following case: a tail of unit chord, tail length I of 3.2 
units, and Mach number 1.7. The result is shown on figure 
39, together with a sketch of the boundary conditions corre- 
sponding to the particular case under study. 



Figure 39. — Response in lift of two-dimensional tail to normal velocity 
of two-dimensional vortex system. Af„=1.7. 


It is to be noted first in figure 39 that the section lift becomes 
infinite when the leading edge of the tail reaches the position 
of the starting vortex (t=B in fig. 39). This is not unex- 
pected, since at the instant before the leading edge penetrates 
the vortex there exists on the tail an infinity in upwash; 
whereas just after the penetration the effect of the infinity 
in up wash is just canceled by an opposite effect caused by 
the infinity in downwash. It is interesting to observe, how- 
ever, that a like result does not occur as the tail trailing edge 
leaves the position of the starting vortex (t—D). This can 
be explained by noting that the sound wave emitted when 
the leading edge penetrates the starting vortex serves to 
broadcast the existence of the infinite discontinuity in normal 
velocity over a finite region of the tail. The distribution of 
the effect of the discontinuity thereby modifies its influence 
at any one point. Next, observe that the lift changes sign 
at O, which is very closely the time when the tail chord is 
situated symmetrically with respect to the starting vortex. 
As can be seen on the flgure, the tail chord experiences very 
nearly equal amounts of upwash and downwash at this time. 
Finally, note that at t=E the lift on the tail disappears as 
the tail moves completely beyond the region of disturbances 
created by the moving vortex at f==0. This fact is consistent 
with the statement made earlier that the final value of the 
lift on the tail should be that due to the steady-state down- 
wash field of the forward wing, for the downwash behind a 
moving two-dimensional vortex in steady supersonic flow is 
zero. It is instructive to remember that the same fact is 
true for the lifting two-dimensional wing in a steady super- 
sonic flow, so that no inconsistency arises from this source by 
a superposition of solutions of the type shown in figure 39 as 
a closer approximation of the response of the tail to a step 
change in angle of attack of the forward wing. The develop- 
ment of such a procedure will be outlined in a later section of 
this paper. 

APPLICATION OF METHODS OF GUST ANALYSIS 

It will no doubt have been noted by the reader that the 
subject of inquiry in the preceding section, namely, the prob- 
lem of calculating the growth in lift at the tail caused by its 
penetration of a region of disturbed flow created by a forward 
lifting wing, bears in many respects a marked resemblance 
to another more familiar problem. That problem is the 
calculation of lift on a wing which penetrates a gust of pre- 
scribed spacewise normal velocity distribution. Indeed, the 
two problems are identical with one exception: in the gust 
problem it is generally assumed that the normal velocity of 
the gust at each station in space remains constant with time ; 
whereas in the wing-tail interference problem this is not the 
case. If it can be shown in the wing-tail problem, however, 
that at each station in space the variation of normal flow 
velocity with time can be safely ignored during the time 
interval required for the tail to pass the station, then the use 
of a gust-type analysis is valid, and an important advantage 
accrues. One can then calculate in a relatively simple man- 
ner the growth of interference lift on any tail surface for 
which the response in lift to a sharp-edged gust of uniform 
intensity is Imown. The so-called gust function has been 
calculated for the two-dimensional case (ref. 3), the rectangu- 
lar plan form (ref. 6), and the triangular plan form having 
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supersonic edges (ref. 5) . Thus, although it is still necessary 
to retain the two-dimensional nature of the downwash field 
created by the forward wing, the gust analysis, if valid, pro- 
vides a means whereby two important classes of three- 
dimensional tad surfaces can be treated. The object of this 
section will therefore be to investigate the validity of the use 
of a gust-type analysis in place of an exact analysis. 

The extent of the approximation involved in the gust-type 
analysis can be best visualized in the two-dimensional case. 
Further, since the exact solution for this case has been ob- 
tained in the preceding section, a comparison of the two re- 
sults provides a means of assessing the magnitude of the 
error introduced by the approximation to the flow. Compare 
first, then, the boundary conditions at the two-dimensional 
tail for the exact and approximate cases. On figure 40 
these conditions are shown schematically for the Mach num- 
ber, tail chord c, and tail length I which correspond to the 
numerical results given in figure 39 of the preceding section. 
For the exact case, the normal velocity at the tail is seen to 
vary inversely as oc or < along rays from the original position 


X 



Figure 40. — Comparison of exact and approximate normal velocity 
distributions at the tail location. 


of the moving vortex. For the approximate case, we fix the 
correct variation of normal flow velocity at the tail leading 
edge and require that at each station x the normal velocity 
existing at the leading edge remains the same for all time 
thereafter. It then is clear from inspection of figure 40 that 
at least for the particular conditions chosen for study, the 
differences in normal velocity experienced by the tail for the 
two cases are not excessive. 

By virtue of the approximation, which fixes with respect 
to time the normal flow velocity at each station in space, the 
problem of calculating the lift at the tail is reduced simply 
to a superposition procedure entirely analogous to that 
already described in the first part of this report for the case 
of a wing performing arbitrary maneuvers. Here, however, 
rather than responding to a series of step changes in angle of 
attack or pitching velocity, the tail responds to a series of 
sharp-edged gusts, so that the Duhamel integral, expressing 
lift at the tail as a function of time is written as. 


C 

In equation (40) the quantity ^ (t") is the gust function 

for the two-dimensional surface, given in reference 3, and w{t) 
is the variation of normal flow velocity experienced by the 
tail leading edge. Substituting x=l—Mot in equation (35), 


w(t)=Mo 


A4> 

2 tt I {I — Mot) 


(41) 


With the use of equations (40) and (41) and the results of 
reference 3 the interference lift at the tail may be computed 
in a straightforward manner. As in the exact problem, how- 
ever, integration difficulties prevent the presentation of a 
solution in closed form in all but the first time interval 
I l-hc 

, - <t< The solution in quadrature form is given 

Mo+1 Mo+l ^ ^ 

in Appendix B. These results have been utilized to compute 
the variation of tail lift with time for the same set of condi- 
tions chosen in the preceding section, namely; Mach number 
1.7, tail chord of unit length, and tail length of 3.2 units. 
Results of these computations are compared with the exact 
results in figure 4 1 . 



f 


Figure 41. — Comparison of exact and approximate tail lift responses. 

M.= 1.7. 
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It is readily apparent from inspection of figure 41 that at 
least for the conditions chosen, the gust-type analysis pro- 
vides an excellent approximation of the variation of tail lift 

I l-\-c 

with time over the entire time interval, , <t< • 

’ Mo-fl M„—l 

We may therefore proceed to the three-dimensional cases 
with reasonable expectation that the solutions resulting from 
the gust approximation will be of acceptable accuracy — 
again, at least for sets of boundary conditions corresponding 
to those used above. 

The analysis for the rectangular and triangular tail plan 
forms proceeds exactly as that for the two-dimensional case 
just discussed. For the rectangular plan form, in fact, the 
two-dimensional solution constitutes a major part of the 
total lift, and it is only necessary to add to the results given 
above the contribution of the tail tip regions. The gust 
function for the tip regions of the rectangular plan form is 
presented in reference 6 and is of such simple form that the 
integration of equation (40) may be done quite easily. Not 
so for the triangular plan form (ref. 5) ; for this case graphical 
or numerical techniques are again necessary to obtain results 
for specific cases. Solutions for both the tip regions of the 
rectangular tail and for the triangular tail are given in 
Appendix B. Numerical results, illustrating the nature of 
the solutions, are given in figure 42. The boundary condi- 
tions for these results have been chosen to be identical to 
those picked for the two-dimensional results of figures 39 
and 41, so that the results of figure 42 for the rectangular 
tip regions may be directly added to those of figure 41 to 
give the total lift for the rectangular tail. Note in figure 42 
that aspect ratio A appears simply as a multiplying factor, 
so that these results are applicable for tails of any aspect 
ratio, subject only to the limiting conditions prescribed in 
the figure. The condition on the rectangular tail insures 
that the Mach line from one edge does not intersect the 
opposite edge; that on the triangular tail insures that the 
leading edge is supersonic. Note also in figure 42 that unlike 
the case of the two-dimensional tail, no infinity appears in 
either the lift contribution of the rectangular tip regions or 



/ 

Figure 42. — Approximate response in lift of wide triangular tail and 
tip regions of rectangular tail to normal velocity of two-dimensional 
vortex system. Mo=1.7. 


the lift of the triangular tail when the tail first penetrates 
the position of the starting vortex. This, of course, is attrib- 
utable to the fact that the first penetration in either case 
is made by a surface of infinitesimal span. 

APPLICATION OF STRIP THEORY 

In the preceding section it was shown that approximate 
response functions could be obtained for those three-dimen- 
sional tail surfaces for which the indicial gust function was 
known. There is stiU another method by which the responses 
of certain other classes of three-dimensional tail surfaces 
may be derived, and that is by the use of strip theory." 
It has been shown by a number of authors (refs. 4, 25, and 
26) that the use of strip theory is exact in calculating the 
lift, pitching moment, and rolling moment of surfaces having 
all supersonic edges and straight trailing edges normal to the 
free-stream direction. Further, no restriction need be placed 
on the generality of the normal-flow-velocity distribution 
encountered by these plan forms. Thus, for example, in the 
present connection, the lift on the wide triangular tail in 
response to either the two-dimensional or three-dimensional 
vortex field can be computed exactly if the responses of the 
two-dimensional tail to these fields are Itnown. 

When restrictions are placed on the nature of the normal- 
flow-velocity field, still other types of plan forms may be 
treated without approximation by the use of strip theory. 
An obvious example is the response of the apex-rearward 
wide triangular tail to the two-dimensional vortex field. 
Clearly, since the loading at no point on the tail can be 
influenced by disturbances from the edges of the tail, the 
integrated lift at every section is just the lift on a two- 
dimensional tail having the same tail length and chord as 
the section. The same would be true for any tail having 
supersonic edges and a straight leading edge normal to the 
stream direction, so long as the normal flow velocity is 
invariant with the spanwise direction. 

As simple examples of the use of strip theory in the present 
problem, formulas will be developed from which the exact 
responses in lift of apex-forward and apex-rearward triangular 
tails to the two-dimensional vortex field may be obtained. 
Consider first, however, the more general plan form pictured 
in figure 43, in which the dimension I is defined as the distance 
between the moving vortex and the foremost point of the 
tail, and the xy coordinate system is fixed in space. 

The use of strip theory (which, it should be cautioned, is 
not exact for the plan form shown in fig. 43) enables one to 
formulate immediately the lift on the finite-span tail in 
terms of the two-dimensional response as follows: 



o \bto/2 


(42) 


where 



{x,t ; y) denotes the loading at a point on a two- 


dimensional tail which is subjected to the same normal flow 
velocity variation as a section c{y) of the three-dimensional 
tail. The inner integral in equation (42) is recognized, 
however, as being just the two-dimensional section lift 


» The term “strip theory” is used here in its usual sense; that is, the integrated lift at a 
spanwise station of the three-dimensional surface is assumed to be the same as that on a 
two-dimensional surface having the same chord length and undergoing the same motion. 
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Figure 43. — Definition of notation used in strip theory analysis. 

coefficient multiplied by the local chord c{y), so that equa- 
tion (42) becomes, 


c{y)Ci{l-]-i{y), c{y),t)dy (43) 


where the quantities l-\-^{y), c(y) replace I and c, respectively, 
in the two-dimensional-tail lift response previously derived. 
The formulas for the apex-forward and apex-rearward tri- 
angular tails are simply special cases of equation (43). For 

2c 

the former case, c{y) is given by y, and after a change in 
variables, y—^ V, equation (43) becomes, 



VCi(l+Co{l-y), 


V^oy f) d?7 


(44) 


For the apex-rearward triangular tail, k{y) is zero, so that, 


a 




-nCo, t) dv 


(45) 


Equations (44) and (45) have been utilized, together with 
the results for the two-dimensional-tail lift response to 
compute the responses in lift of apex-forward and -rearward 
triangular tails to the two-dimensional vortex field. Results 
in integral form are given in Appendix C, and numerical 
results for tails of unit chord, tail length 3.2 units, at Mach 
number 1.7, are shown on figure 44. 

The result for the apex-forward triangular tail is compared 
with the approximate result derived for this wing in the 
previous section by use of a gust-type analysis. Again, as 
in the two-dimensional case, it is evident that the gust-type 
analysis provides an accurate estimate of the tail lift response 


over the entire time range. 


<t< 

M„+l - - Mo- 


0 

l' 


SUPERPOSITION OF ELEMENTARY SOLUTIONS 



Figure 44. — Exact response in lift of apex-forward and apex-rearward 
wide triangular tail to normal velocity of two-dimensional vortex 
system. M„=1.7. 


velocity field created by a forward wing of finite chord. In 
this section a procedure for accomplishing this end is de- 
veloped and the result of the analysis applied to a special 
case. For the sake of simplicity, the analysis is again 
restricted to one dealing with two-dimensional forward 
surfaces, so that the normal-velocity field encountered by 
the tail is uniform in the spanwise direction. However, no 
restriction is placed on the span of the tail. 

In the previous sections solutions have been obtained for the 
lift at the tail in response to a single vortex pair. The 
boundary conditions for the vortex pair are exhibited in the 
plane as a triangular sector having a imiform potential 
difference A<i> within the sector (see fig. 36). We wish now to 
distribute these sectors over the region traced out in the xt 
plane by the forward wing and its wake in such a way as to 
build up the variation of A4> in these regions corresponding to 
prescribed normal-velocity patterns over the wing. The 
sum of the responses in lift of the tail to each of these slabs of 
A4> is then the response to the normal-velocity field created 
by the wing. It is convenient to treat the effects of the wing 
and wake separately, and we consider the wing first. Place 
a vortex pair at an arbitrarj^ chordwise position of the wing 
with origin at f=0. The response in lift at the tail to the 
single vortex pair is a function of time t and a characteristic 
distance, chosen to be the distance n between the moving 
vortex and the tail leading edge. This of course is the solu- 
tion derived previously with I replaced by Now place 
another vortex pair of opposite strength with origin at the 
same position in space as the first pair, but shifted an incre- 
mental tim,e Ari later. As shown in figure 45 there remains 
a strip of chord MoAtj having unit A4> within the strip. 

The response in lift at the tail to the remaining strip at a 
fixed time t is, 


It has already been noted in a previous section that the 
response in tail lift to the normal-velocity field created by a 
single vortex loop could be viewed as a fundamental result 
in this sense: The solution can, in principle, be used as the 
elementary function in a superposition procedure designed 
to determine the response in lift of the tail to the normal 


AOl(ix, f) = Cr,(g, t) — CL{ix—M oAy, t—Ari) (46) 

Letting the changes in and t approach zero, equation (46) 
may be cast in the form of the directional derivative, 

^ (m, 0 =M„ ^ (m, 0 + ^ (m, t) (47) 
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Figure 45. — Boundary conditions for response in lift of tail to an 
elementary strip of unit A4>. 


This is the response in lift to a strip of unit A$ situated a 
distance p from the tail. The response to any other strip 
starting at a later time but at the same distance n from the 
tail will of course be identical to this result. We may there- 
fore build up the response to an arbitrary variation of A# in 
the strip simply by use of the Duhamel integral ; that is, 

(48) 

Let t — ti = T, and perform the indicated differentiation, noting 
from physical considerations that A4 >(m, 0) = C'l(m, 0) = 0. 

^ (m, 0=J^ \_Mo ^ (y, r) + ^ (y, r)] A$ (m, t-r)dr 

(49) 


The influence of the wing on the lift at the tail is then ob- 
tained simply by a summation of incremental strips as p 
goes from I to l-\-c„^. Note that the increment Ap is Mc^rf. 

a4> (p, t — r) dr 

(50) 



Now consider the influence of the wake. The potential dif- 
ference A# in the wake has the characteristic that at each 
station in space the A$ existing at the trailing edge of the 
wing remains the same for all time thereafter. The distribu- 
tion of A^ in the wake may therefore be built up by elemen- 
tary vortex pairs whose origins are placed at the trailing edge 
of the wing. Since these elementary sectors are all the same 
distance I from the tail, the response in lift at the tail to each 
of the sectors is of the same form, and the sum of these re- 
sponses at a fixed time t may again be found by use of the 
Duhamel integral, 

GL{f)wAKB—-^^ ^ Ct{l, t)A4>(Z, t r)<fr 



r)Ai(f, t — t) dr 


The sum of equations (50) and (51) is then the desired re^ 
sponse in lift at the tail to a prescribed A4> distribution on the 
forward wing. 

Ci(f)= Jj dnj^ in, (m, t)Ja#(m, f—T)dT-f 

^‘^CS,r)mi,t-r)dr (52) 

We consider next a special case for which the form of equation 
(52) may be simplified considerably. For this case the 
loading on the forward wing is prescribed to be uniform 
in the chordwise direction and to remain constant with 
time. The specified boundary condition may be interpreted 
either as the exact case for a forward wing which starts to 
deform at t=0 in such a way as to give rise to a uniform 
loading over its surface, or as an approximation to the res- 
ponse of a flat forward wing to a step change in normal 
velocity over its surface. It is the latter case which is of 
real interest here, and of course the approximation introduces 
some error into the variation of A$ on the wing and in the 
wake. However, as can be seen on figure 46, a comparison 
of the exact and approximate A«l> distributions shows that 
the differences between the two are not great for the case 
chosen, and would be even smaller for higher Mach numbers. 
It is expected that the appro.ximation should yield acceptable 
results for Mach numbers greater than about 1.4. 



Figure 46. — Comparison of distributions of for two-dimensional 
wing and wake due to uniform loading and uniform normal velocity 
on wing. 

The prescription of uniform loading on the forward wing is 
seen on figure 46 to amount to the following conditions on 
A$ on the wing: 


where 


A$(/i,f)=H 0<t<t 

=kt t>t 


_ l+co^—n 
Mo 


(53) 


(51) 


It follows that the condition on A$ along the trading edge 
of the wing must be, 
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A^(l,t)=kt 




-k^ 

M„ 


Co 

t> irf 
-Mo 


(54) 


Performing the required differentation of A$ with respect to 
{t—r) and substituting in equation (52), we get 

7i h 

k ^l_c^GL{l,r)dr (55) 


M, 


The inner integral h in equation (55) may be evaluated 
immediately and gives 




(56) 


For the inner integral h, we make use of the definition for 
the differentiation of a definite integral, that is, 

f7i(M,r)dr-(7i(M,f)0+ 

(57) 

Substituting equations (56) and (57) in equation (55) we 
have, 

^<7i(M,f)]+fcj|^^(?i((,r)dr (58) 

Mo 

Let 

f .CLitl,T)dT = F{fJ.,f) 

Jt-t 

Then in equation (58) 

Ci+c.„ A 

r Oi(f+c„^,r)dT— r c^Gi^{l,T)dT 
Substituting equation (59) in (58) there remains only, 
6i(t)=-^J^ CL(fj,,t)dfi 


(59) 


(60) 


This result may be easily verified by interpreting the integral 
in equation (60) physically. The integral specifies that we 
sum the effects of elementary sectors, all of the same strength 
in A$, and whose origins are at f=0 within the boundaries of 
the wing leading and trailing edges. A graphical summation 
of several such sectors will show they do indeed build up the 
A$ distribution over the wing and wake prescribed in figure 
46. 

Equation (60) has been utilized along with the results of 
the previous sections for for the two-dimensional tail 

and the rectangular tail tip regions to approximate the 
response of the tail to a step change in angle of attack of a 


forward surface of finite chord. Results of computations 
for a wing chord of 2 units, tail chord of 1 unit, tail length 
f of 3.2 units, and Mach number 1.7 are shown in figure 47. 



Figure 47. — Approximate response in lift of two-dimensional tail and 
tip regions of rectangular tail to step change in angle of attack of 
two-dimensional wing. M„=1.7. 


Note that the main difference for the two-dimensional result 
from the result for the response to a single vortex pair is the 
disappearance of the infinity in lift when the tail penetrates 
the position of the first starting vortex. Reference to figure 
46 will show that this is attributable to the fact that the 
first vortex penetrated by the tail is of infinitesimal strength . 
Again, as in the single-vortex example, the response in lift 
for a finite-span tail satisfying the restriction >1 may be 
obtained from the results of figure 47 by direct addition of 
the two variations shown in the figure. 

It should perhaps be noted that the form of equation (52) 
permits of still other simplifications which should yield 
acceptable results for special cases. For instance, when the 
chord of the wing is small compared to the tail length, it may 
be expected that the influence on the build-up in lift at the 
tail of disturbances on the wing should be minor compared 
to the influence of disturbances in the wake. It is seen in 
equation (52) that this condition is equivalent simply to 
ignoring the first integral in comparison to the second. The 
simplification, of course, represents a considerable reduction 
in computational labor, and, in addition, the form of the 
remaining integral is simple enough so that the exact varia- 
tion of A<i> at the trailing edge of the wing corresponding to 
the prescribed normal-velocity distribution on the wing may 
be used with relative ease. For example, for the case of 
uniform normal velocity on the wing, the variation of A$(f,<) 
at the wing trailing edge may be derived from the results of 
reference 3, and is given by 


AHl,t)=Voa, 


Voa 


cos ' 



^0^ </.< 

Mo — 1 




(61) 


= 0 , 


Mo-l~ 
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Upon substitution of these expressions in equation (51) and 
performance of the required integration, the solution yields 
in any case the exact response of the tail lift to the influence 
of the wake, and, in cases where the tail length is large 
compared to the wing chord, a reasonable approximation ol 
the total response in lift of the tail. The response in lift 
of the tail to, say, indicial pitching motion of the wing, or 
for that matter, to any indicial, motion of the wing for which 
the variation A$(Z,i) is known, can be approximated by 
the same technique. 

CONSIDERATION OF THREE-DIMENSIONAL FLOW EFFECTS 

In all of the analysis presented thus far, attention has 
been focused solely on problems in which the flow field 
caused by the forward wing is two-dimensional. This 
effectively limits the results to cases wherein the plan form 
of. the forward wing is essentially rectangular, with span 
sufficiently large so that the flow from the tips of the wing 
does not significantly affect conditions in the vicinity of the 
tail. In order to treat more general cases, including cases 
wherein the plan foi’ms are of nearly equal span, it becomes 
necessary to consider in detail the effects on the tail of the 
remainder of the vortex loop — the trailing vortices. 
Inevitably, the admittance of spanwise variations of normal 
velocity into the flow-field picture complicates the analysis 
considerably, and it appears that exact solutions for the 
tail lift response will be exceedingly difficult to obtain. 
However, the change in the physical situation caused by the 
addition of the trailing vortices can be fully described, and 
this description will be seen to suggest simplifying approxi- 
mations which serve, in turn, to make the problem of cal- 
culating the tail lift response more tractable. 

THE PLOW FIELD 

Consider first the flow-field caused by the motion of the 
complete vortex loop. The boundaiy condition character- 
istic of the vortex loop is that within the confines of the 
moving, starting, and trailing vortices, the jump in potential 
A$ is a constant across the plane containing the vortices, 
and zero elsewhere in the plane. Thus it is clear that the 
two-dimensional vortex system considered earlier can be 
converted to the three-dimensional system by subtracting 
off regions of constant A$ of semi-infinite spanwise extent, 
as shown in figure 48. 
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Figure 48. — Superposition procedure for converting vortex system 
from two- to three-dimensions. 


Obviously, by reason of symmetry, we need consider the 
effect of only one of the added trailing-vortex loops, and, 
further, it is sufficient to consider the effect only within the 
shaded region of figure 49. 



Figure 49. — Boundary conditions in the 2 = 0 plane for a trailing 

vortex loop. 

In figure 49 the circular region, denoted region I, is the 
domain in the plane of the vortices within the sound wave 
emitted al <=0 by the motion of the trailing vortex. This 
region extends from values of x betiveen -\-t and — f, and 
constitutes the interval within which the flow is both un- 
steady and three-dimensional. For values of x<^ — t, the 
flow has reached its steady-state distribution, and is therefore 
stationary with respect to an observer traveling at the speed 

of the moving vortex. The lines y= ^ setting off 

this steady-flow regime are of course the Mach lines 
emanating from the juncture of the moving and trailing 
vortices. All points to the left of the shaded region are 
clearly unaffected by disturbances from the trailing vortex 
itself, and hence the flow there is just the mirror image of the 
two-dimensional flow considered earlier. Points to the right 
of the shaded region are completely unaffected by the action 
of the trailing-vortex loop, and there the fluid is undisturbed 
by the addition of the loop. 

The problem of calculating the potential and velocity 
field associated with the trailing-vortex loop has already been 
considered for the subsonic speed range by Lomax, Heaslet, 
and Fuller in reference 24. It is a relatively simple matter to 
extend these results to the supersonic range, and since the 
procedure for doing so is plainly developed in reference 24, 
only the final results that are of interest here will be given. 
Thus, for the regions marked I, II, III in figure 49, the 
normal velocity w induced in the plane of the vortices by the 
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addition of the trailing vortex loop is given by : 


HKSPONSE IN LIFT OF TAIL TO NORMAL VELOCITY OF TRAILING- VORTEX 

LOOP 


- n n\ r - ^{x+Moty—fiY 
y — , _x{x-\-Mot) y{x-\-M„i) 




xy 




2x y{x+Mot) 


1 /< 0 , ^W-^2%x{x-^Mpt)'^A.T\_ x(x+Mpt) 

^{x+Motf—0^y‘^ , 


y{x-\-Mot) 


xy 


-w(II) = + 


A-l» Mp-^Ji^—x^ ■\j{x+Mptf—fi'‘y^ 
2-ir x{x-\-Mpt) 2ir y{x+Mpt) 


h (62) 


A$ ^{x+Mpty—0^y^ 

-W{111)= ~ y(^^+Mpt) 


Several points of interest may be deduced from inspection of 
equation (62). First note that, as it should, the solution 
reduces to the t-wo-dimensional normal velocity with just a 

change in sign along the left-hand boundary (^= — -yjt^—7?, 
and to zero along the right-hand boundary (y= + -\jf—x‘, 

='=+C^> w. >->-«.')■ 

Second, it appears that the contribution to the normal 
velocity made by the trailing-vortex loop may be viewed as 
being composed of two distinct parts; one, a two-dimensional 
contribution, equal but opposite in sign to the normal velocity 
contribution of the original two-dimensional system and 
existing everywhere within the region y<0, —t<Cx<Ct, and 
second, a three-dimensional contribution that is antisym- 
metric with respect to the axis y=Q and that exists only 
within the confines of the starting sound wave and Mach line 
traces, that is, within the shaded region of figure 49. It 
should be noted that the antisymmetric contribution is zero 
along both the left-hand and right-hand boundaries, and 
contains no singularities other than at y=0, x<C.O, the path 
of the trailing vortex itself. 

In order to obtain the normal-velocity distribution due to 
the complete vortex loop, it is now only necessary to combine 
the results given by equation (62) with the two-dimensional 
contribution (eq. (35)), and with the contribution of the 
trailing-vortex loop originating at y=s (see fig. 48). How- 
ever, we are interested in the effect of this complete loop on 
the lift at the tail and, inasmuch as the response to the two- 
dimensional vortex system has already been investigated in 
previous sections, it is more to the point simply to continue 
by investigating separately the effect on the tail of only one 
trailing-vortex loop. The resulting solution for tail lift re- 
sponse can then be doubled to account for the influence of 
both trailing-vortex loops, and directly added to the result 
found previously for the response to the two-dimensional 
vortex system. 


In view of the complicated nature of the flow caused by 
the trailing-vortex loop, it appears that analyses leading to 
exact solutions for general tail plan forms will be confronted 
by formidable diflflculties. Thus, for practical purposes, 
simplifying approximations will no doubt become necessary. 
There is one type of plan form, however, which can be 
treated with a minimum of approximation and that is a 
class of supersonic-edged wings. It has already been men- 
tioned in a previous section that for the supersonic-edged 
plan form having a straight trailing edge normal to the free- 
stream direction, the use of strip theory is exact in calcula- 
tions for the lift due to arbitrary deformations of the surface. 
Hence, for this type of plan form one can calculate the lift 
at a given spanwise station as though the station were a 
two-dimensional surface undergoing the same motion. The 
reduction of the problem to one involving only two spatial 
dimensions then opens the way for use of the three-dimen- 
sional steady-flow analogy to solve the problem in the same 
manner as was done previously for the response of the 
two-dimensional tail to the two-dimensional vortex system. 

As an illustration of the use of these principles, we con- 
sider two cases involving the response of a tail of triangular 
plan form; first to a vortex system having a span greater 
than, and second a span less than the span of the tail. 

Case 1 ; 2mc„<s. — 

Referring to the notation of figure 50, let us fix attention 
to a spanwise station y, and examine the normal velocity 
distribution encountered by chord Cj (y) of the tail as it passes 
through the region of disturbed flow caused by the left-hand 
trailing-vortex loop. For the particular chord chosen, the 
geometry of its experience in time and space is illustrated 
in figure 5 1 . 

As seen on tlie figure, in virtue of fixing y, the boundaries 
of the sound wave appear as one branch of a hyperbola. 
Within the hyperbola, the flow is of the type defined by 
w(I) of equation (62) ; whereas to the left of the hyperbola 
but within the trace x=Py~Mpt, the flow is constant along 
lines parallel to the leading edge of the section, being of 
type III of equation (62). Elsewhere, the flow is zero. 
Now, in view of the validity of the strip-theory formulation 
we consider the chord Ci{y) as being the chord of a two- 
dimensional tail, subjected to the normal-velocity distribu- 
tion w{x, t; y), where y is treated as a constant. Then, 
invoking the three-dimensional steady-state analogy, the 
problem is to find the steady-state section lift as a function 
of t (the coordinate t is now the stream wise axis) on the 
equivalent sweptback tail shown in figure 51, which is sub- 
jected within the shaded boundaries to a normal-velocity 
distribution given by equation (62). 

Up to this point, the formulation of the problem has 
involved no approximations. Now, however, rather than 
attempt to solve the above problem e.xactly (which is in 
truth even now a formidable exercise in integration) we 
introduce a simplifying approximation based on the nature 
of the normal-velocity distribution. It has been mentioned 
that the antisymmetric part of the normal velocity caused 
by the trailing-vortex loop contains no singularities other 



USE OF INDICIAL FUNCTION CONCEPT IN ANALYSIS OF UNSTEADY MOTIONS OF WINGS AND WING-TAIL COMBINATIONS 31 



response to normal velocity of trailing vortex loop. 



Figure 51. — Boundary conditions at a spanwise section of the tail. 

than those on the path of the trailing vortex itself. It can 
also be said that with respect to an observer moving with the 
chord Ci{y) the magnitude of the antisymmetric flow grows 
uniformly from zero at one boundary of the sound wave 
to the steady-state value at the other boundary. These 
statements are illustrated in figure 52 for the region in x, 
t, occupied by the tail chord Ci{y) of figure 51. 

Thus, from inspection of figure 52 it can also be stated that 
the response in lift Cj(^+?i(2/), Ci(y), t) must start from zero 
at time ti{y) (see fig. 51), and also grow uniformly to its 
steady-state value, which is attained at time tziy). More- 
over, it can be shown that if w is continuous along lines 



Figure 52. — Normal velocity distribution duo to trailing vortex loop 
at a spanwise station of the tail. 


parallel to the leading edge the slope of the response curve 
must be zero at both times h and < 2 - 

The steady state lift 'ci{l+^^{y), Ci{y)) which is attained at 
time <2 can be computed exactly, and the times ti and U are 
obtained easily from the geometry of figure (51). It there- 
fore seems reasonable simply to approximate the response 
curve by a variation that grows uniformly with time and 
fits the end conditions correctly in magnitude and slope. 
One such curve can be written in the form of a cubic equation 
in t, and gives, 

where, for the conditions of figure (51), 

^yy_ Mo[i+uy)]~^[i+mr-^Y 


,,, Mo(i+co)+.ja+cor-^y 

hW — ^2 

Ul+M 

l*y (if7.)-VI'+£.W)‘-n¥+% 74^} 

The form of equation (63) can be viewed as being typical 
of the response of any section c(y) of the tail. However, as y 
takes on all values from one tip to the other, the equations for 
ti, ti, and Cl may change from type to type in different span- 
wise intervals, depending on the Mach number and con- 
figmation geometry under consideration. Referring, for 
example, to figure 50, it is clear that a section of the tail near 
the right-hand tip may be completely untouched by the 
normal velocity caused by the trailing vortex; whereas sec- 
tions closer inhoard may be touched on only a portion of 
their length. All told, there are five such possible variations 
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in the boundary conditions and these are listed in Appendix 
D, together with the correspondiug values of ti ,tn, and Cj. 

Having determined an expression for the response in lift 
at an arbitrary section of the tail, we may now integrate 
and fin d the combined response of all sections of the tail to 
the given normal- velocity distribution. Thus, in a manner 
similar to that described previously in connection with cal- 
culations for the response to the two-dimensional vortex 
pattern, we may write, for the response to one trailing vortex, 

1 f rmc„+R 

Cl {y) Cl (Z+^I (y), Cl (y), t) dy+ 

J ^2mc A-R 

C2{y)ci{l+k2{y),C2(y),t)dy I (64) 

where the symbols in equation (64) have been defined in 
figure (50) . In order to illustrate the nature of the result for 
the response in lift to the trailing-vortex loop, calculations 
based on equations (63) and (64) have been carried out for 
the triangular tail shown on the side of figure 53. The 
results account for the effects of both exterior trailing-vortex 
loops and may be directly added to the response to the two- 
dimensional vortex system to give the response to the com- 
plete loop shown in figure 53. The separate responses and 
their combination are shown for comparison on figure 53. 


where the variation Ci{l-\-^{y), c(y),t) is again given by 
equation (63). 



Figure 54. — Definition of notation used in calculating section lift 
response to normal velocit.v of trailing vorte.\ loop. 



Figure 53. — Response in lift of wide triangular tail of aspect ratio 
4 to normal velocity of complete vortex loop. Mo=1.7. 


Case 2: 2mCo^s . — 

When the span of the trailing vortices is smaller than that 
of the tail, the development outlined above must be amended 
slightly. The effect of the antisymmetric part of the normal 
velocity-distribution may be accounted for, however, by the 
same technique as described previousljq and referring to 
figure 54 for definition of the notation, the response in lift of 
the tail to this contribution may be written immediately as, 

J ^ ci{y)ci{l+^i(y),ci{y),t)dy+ 

C2(y)ci(l+Aj),c2(y),t)dy V (65) 

Jmc^—R J 


Now, it is also necessary to consider the influence of the two- 
dimensional contribution to the normal velocity which exists 
in the region y<0, — <<x< + f. But note, the response in 
lift of the tail to this contribution is simply equivalent to 
minus one half the response of a triangular tail of chord 
R R 

— and tail length l+Co to the two-dimensional vortex 

m ° m 

pattern, and, inasmuch as the exact result for this case has 
been obtained in a previous section, it is sufficient to indicate 
the result in coefficient form as 

where the subscript 2 indicates that the response to the two- 
dimensional vortex pattern is to be taken. The sum of 
equations (65) and (66) then gives the total response of the 
tail to one trailing-vortex loop. To obtain the response to 
the complete vortex loop, this result is then doubled and 
added to the response to the two-dimensional pattern in the 
same manner as was done for case 1 . 

SUPERPOSITION OF ELEMENTARY SOLUTIONS 

Finally, having obtained a solution for the lift at the tail 
induced by the normal-velocity distribution of the complete 
vortex loop, we wish now to use the solution to build up the 
response to the downwash caused by a three-dimensional 
forward wing. The analysis for accomplishing this task 
closely parallels that described previously for the case of the 
two-dimensional forward wing, and therefore only the final 
result will be given here. Thus, denoting the elementary 
solution by CA,s,t), where m is the distance from a moving 
vortex to the tail leading edge and s is the semispan of the 
vortex loop, we get, (see figure (55) for definition of notation) 
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fi(() ri 

A^(M,s,t—T)dr-i-J dsj 


dfids 

ds 


(l,s,T)AMl,s,i-T)dT (67) 



of tail to prescribed Alt’ distribtition on three-dimensional forward 
wing. 

As in the two-dimensional case, the fii’st term in equation (67) 
accounts for the variation of A<i> on the wing, the second, 
for the variation in the wake. It is to be noted that in the 
development leading to equation (67) it is assumed that the 
trailing edge of the forward wing is straight and normal to 
the free-stream direction. However, wings having other 
than straight trailing edges can be treated by the same 
equation by considering as part of the wing points in the 
wake between the wing trailing edge and the rearmost 
spanwisc axis passing through the wing. 

APPLICATION OF REVERSE FLOW THEOREM 


tion, and boundary conditions for the flow over the surfaces 
are specified only for <>0, being zero for <<0. After the 
combination has moved for a time T in both forward and 
reverse dflections, the relationship of the wings and their 
respective coordinate systems are as shown in figure 56. 



Ficcre 56. — (Coordinate systems in forward and reverse motion. 


In figure 56, the subscript 1 denotes the axis system in for- 
ward motion, the subscript 2, the axis system in reverse 
motion, and the wing, wake, and tail have been given the 
symbols a, b, and c, respectively. As may be verified 
from examination of the figure, the x and t axes are related 
by the equalities, 


X2=—Xi+'i—MoT 
ti= — ty-\-T 


(68) 


and wc furth.er specify that 1 / 2 =— 2/i- The general reversal 
theorem may then be written in the form (see again, ref. 27), 




(I.) 


[Pa,(x^,Vl,U)'Wa.^(Xi,yi,t\) + 


We have devoted several sections to the derivation of 
methods whereby the tail lift response can be calculated for 
certain restricted cases. The possibility still exists, however, 
that these results can be extended to include still other cases 
which do not lend themselves as readily to analysis. The 
question may be asked, for instance, whether a relationship 
exists between the response in lift at the tail to the essentially 
two-dimensional normal-velocity field shed by a large-span 
rectangular wing and the analytically more difficult reverse 
situation; that is, the effect of a small-span forward wing of 
fairly general plan form (for which the downwash field is far 
from being two-dimensional) on a large-span rectangular 
tail. Such a relationship can be sought for the more general 
case, involving wing-tail configurations of arbitrary size and 
shape, by application of the reverse-flow theorems for un- 
steady motion developed in reference 27. It will be the pur- 
pose of this section to investigate this possibility. 

To begin, we adopt the notation used in reference 27 and 
also specify that the motion of the wing-tail combination be 
indicial; that is, the combination is made to start from rest 
at f=0 with uniform velocity Vo in forward and reverse mo- 


Vb, (J^i, ^1, h)Wi.,(^i, 1/1, h)+Pc, i:xi,yi,k)Wo, fi)] dxidyi= 

f [Pa,{X2.y2,t2)Wa^(X2,y2,k) + 

Jo 


Pb.^(x 2 ,y 2 , k)Wb,(x 2 ,V 2 , ( 2 ) -Vvc^(x 2 ,y 2 ,h)'^c^(x 2 ,y 2 , ^ 2 )] dx^dyi (69) 


where P is the plan form of the combination in the s=0 
plane. It may be noted the normal-velocity functions, 
Wq, are implicitly dependent on T, the full relationship 

being. 


Wq (X 2 , 2/2, <2) = •Wq ( — X2 -b f— T, — 2/2. — ^2 + T) ' 

Wq(x,, 2/,, U)=w^( — XiV^—MoT,—yu — U+T) - 

i=a, 6, c 


(70) 


Notice further that, since the wake cannot sustain a pressure 
difference, the products and Vb^b^ in equation (69) 

may be set to zero. Then taking a derivative with respect 
to T oi both sides of equation (69), we get. 
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iX ^_^^V<:Muyi,T)Wc^{xi,yi,T)dxidyi + (, > 


i>Wa, 

bU 


^dxidyi+^^ ^^'XC (,, ^^^Va^{^,yi,T)Wa^{x2,y2,T)dXidyi 

XC (y,, 2/2, 1/2,7 ’) ^^XX ^“ 2 (^ 2 ’ 2 / 2 ’ 


+ 


f duff Pcfe, 1 / 2 , ^ 2 ) — M„ 

Jo JJp^^u,) - \ OXi Oti / 


dxidy^ 


(71) 


It is easily verified that a sufficient condition for the disap- 
pearance of the triple integrals in equation (71) is that the 
normal-velocity distributions on the surfaces be restricted to 
those defined by the relations, 


2/i)=/(xi+MJi; yi) 


W 2 {x 2 ,h] yi)=g{x 2 ■^-M„t 2 ] 1 / 2 ) 


(72) 


Another such condition is simply that Wi = 0, or W2=0. 
Assuming henceforth that either or both of these conditions 
are in force, the flow reversal relationship for the wing-tail 
combination becomes, 


J J Pai (xi ,yi , T)Wa^ (xi ,yi , T) dxidyi_ -|- 

JJ Pc-iixi,yi,T)wc^(Xi,yi,T)dxidyi 

Pc,(T) 

-jj Pa2(a-2,i/2, T)wa^(x2,y2, T)dxody2+ 


(73) 


jj Vc^{x2,y2, T)wc^(x2,y2, T)dx2dyo 


A"ow in the present problem, the boundary conditions for the 
flow on the wing and tail are simply that in forward motion 
the normal velocity on the forward wing is uniform while on 
the tail it is zero, and in reverse flow we specify that the same 
situation be true. Obviously, these conditions satisfy the 
restrictions imposed by equation (72), and setting Wa^ and 
equal to unity, equation (73) becomes, 


jj Pci(xi,yi,'ildxidyi= j j pa^{x2,y2,T)dx^dy2 (74) 


Pa,(T) 


APPLICATION OF RESULTS TO DYNAMIC STABILITY 
ANALYSIS 

The pi-eceding sections have been concerned with the 
development of methods whereby the variation with time of 
the lift developed at the tail in reaction to the normal 
velocity shed by a forward wing could be estimated. Now, 
having arrived at a fairly rigorous physical understanding of 
the nature of the variation, we will attempt in this section to 
assess the significance of its contribution to the over-all 
dynamic stability of an aircraft. For simplicity, attention 
is devoted to the single-degree-of-freedom rotary-oscillation 
case. 

It has been noted that the oscillatory lift and moment 
equations developed previously for a single wing are still 
applicable for the wing-tail case with proper interpretation 
of the indicial functions. For supersonic speeds, this re- 
definition was seen to reduce for each indicial function simply 
to the sum of the indicial functions of the wing and tail 
considered as isolated wings and an interference function 
which accounts for the influence of vorticity shed by the 
wing on the development of lift at the tail. Then, denoting 
the interference functions by primed quantities, we conform 
to the notation used in Part I by letting 

a,(r) = njO-^’:(r)' 

C„(r) = r:j/:)-F^(r) 
r’^Jr) = aj/:)-F'(r) I 
r:,(r)=r:jO-F;(r) 


Again introduce the parameters, 




number of half M. A. C. lengths traveled in 
time T 

reduced frequency 


Thus, the response in lift induced on the tail by the 
indicial change in angle of attack of the wing is exactly the 
same as the lift response of the wing due to indicial change in 
angle of attack of the tail in reverse motion. It will be noted 
that the above result applies for wing-tail combinations 
flying at either subsonic or supersonic speeds, and that no 
restriction has been placed on size or shape of the plan forms. 



number of half M. A. C. lengths for inter 
ference function to reach steady-state. 


where the characteristic unit of length c is taken as the wing 
M. A. C. The in-phase and out-of-phase lift and moment 
coefficients for the harmonic rotary-oscillation case then 
follow directly from the form of equations 12 and 14. 
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R.P. j J JO - 


ij' [FiJ^fp)-\-FiJ^<p)-\-F[(‘p)\ sin kipd<p-\- 
^ {F 2 Jjp)-\-F 2 J^<p)^F'^{<p)\ cos, k<pdtp 

* f ° [-^2„(^)+^2j.(v’)+jP’Jip)] sin k<pd(p— 

Jo 


j: 


[Fi^{(p)+Fi^{<p)+F[{<p)] cos k>pd<py 


(76) 


y (77) 


k^ [FiJ^ip)-\rFiJ^ip)-{-F'i{tp)\ sin A;^(Z<p-|- 
F y [F,S‘p)+F,,{<p) + F:{^)\ cos k^pdv 
LP. (yyk{ j + C'„,^(<x„ j + Giy:) - 

k^ [FiJ^(p)-\-Fiy)-\-F'i(ip)\smkipdip 

J' l-^3„,(‘P) + 7^3j'P) + -P’J^)] COS ktpdtp^ J 

In the same way, the reduction of the above results for 
to first-order in frequency follows directly from the 
form of equations (13) and (15). There results, 

[Fiy)-\-Fiy)-\-F[{(p)]dip'^ 

R,P, (yyc^y.y+c^yc.;)+cLs-:) 

/.P.(^)=fc{0.^J.O + ^>»„„(v) + ^;(O- (79) 

J" [F3y)-\-F3y)+F'3{<p)]d<p^ 

With the aid of the first order in frequency results (eq. 
(79)) it is possible to assess the importance of interference 
effects on the damping moment of a wing-tail combi ration. 
Inspection of equation (79) for the out-of-phase pitching 
moment reveals that the interference effects are contained 

in the term and the integral, — J F' 3 ((p)d<p. The 

12 It is presumed in equations (75) to (79) that all of the coefficients have been nondimen- 
sionalized with respect to the same characteristic area and length. 


term is generally relatively small, and may be 

approximated by replacing the steady nonuniform downwash 
at the tail caused by the pitching velocity of the wing by an 
appropriately averaged effective angle of attack at the 
tail. The term C^JtrJ may then be related to Gmy<ra^) 
by the expression, 

(80) 


^IV 

The downwash function -4 has been plotted for certain cases 

2c ... 

in reference 28. For other cases, Ribner indicates in refer- 
ence 29 that a generally suitable expression in terms of the 

V) 

more fully investigated quantity is 


qc aVo 

Thus, combining equations (80) and (81), 


(81) 




(82) 


For pitching in the positive sense the effective angle of 
attack at the tail will in most cases be increased, so that 
equation (82) represents a small but stabilizing damping 
moment. 

Consider now the term — J' FJip)<f‘P- For axes of rota- 
tion well ahead of the tail, the form of the pitching moment 
variation G'„Jip) will vary closely as— ^ G'lJ,<p), wheic It 

may be taken as the distance from the steady-state center 
of pressure of the tail to the axis of rotation. Then, recalling 
the general nature of the curves found for tail lift in previous 
sections, a typical variation G'„ {tp) can be drawn as shown 
in figure (571. 



FiGURi- 57. — Typical response in pitching moment of tail to step 
change in angle of attack of wing. 

The term— J" ° Fj<p)d<p is then seen to be equivalent to 

the net shaded area in figure (57). Plainly, it can be of 
significant magnitude, and since the steady-state ordinate 
GLg{<r'«) tvill most generally be positive, it represents a 
potentially large stabilizing damping moment. 

It is interesting to compare the interference damping 
moment given by other theories in terms of the result given 


IIHI IMailllill ■ I I III! I IHII ■lll■lllll II lllllll 


■■III I nil nil I III nil I III lllllll I III 
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"by the integral J" ° K (<f>)d<p. Most of these theories (c. f., refs. 

12, 28, and 10) are similar in that they recognize that an 
out-of-phase lift at the tail is caused by the fact that a 
finite time is required for the normal velocity shed hy the 
wing to reach the tail. The normal velocity at the tail, 
. which is usually taken to be the normal velocity in steady 
flow, is then usually averaged to form an effective angle of 
attack, and the damping moment cast in a form which may 

he represented in the present notation as — g- 

Now this result may also be interpreted as an area, being 
in fact the area of a rectangle having as ordinate the quan- 

tity Cm (o"aJ. But note that this ordinate is equal to 

aV o “r ~ 

the ordinate Cm„{a'a) in figure 57. Thus, it can be said 
that the classical “lag in downwash factor" is essentially 
correct but depends for its accuracy almost entirely on the 

21 ' 

proper choice of the effective tail length In reference 

12, the length I' is chosen to be the mean distance between 
the tail and the wing. Thus, for rectangular plan forms, 

21' 

and in the notation of figure 34, the effective length ^ is 


■Cir+Cr+2i' 


'aK ““A 


The damping moment is then equal to 
equivalent 


area is com- 


pared with the areaj^ F'z{<p)d(p in figure 58. 



dowmvash factor. 


It can be seen from examination of figure 58 that the simple 
representation of the tail length is in this case surprisingly 

21 ' 

accurate, since the length is so situated that the area 

marked a representing a stabilizing contribution is in large 
part canceled by the destabilizing contribution, area 6. 
Further, it may be argued that this will be true for a large 
number of cases. The reasoning is as follows: Referring to 
figure 58, we see that the use of the simple tail length will 
he essentially correct so long as the number of half-chord 
21 ' . 

lengths — is close to the value of <p at which the pitching- 
c 

moment variation changes sign, for then the areas a and b 


tend to cancel. But the pitching moment will be near 
zero when the tail is so situated in the normal velocity field 
from the wing that it experiences equal amounts of upwash 
and downwash over its surface. Reference to figure 59 
shows that this will occur when the tail is approximately 
midway between the starting vortices from the wing lead- 
ing and trailing edges. 



Figure 59. — Sketch showing tiine at which tail experiences approxi- 
mately equal amounts of upwash and downwash. 


The time t' when this occurs may be calculated to be 



(r+f+l). 


of travel, we get ^ ~ 


and, after converting to half-chords 
^ 2^-hcy-hcr ^^ which is just the value 


of the effective tail length, 

c 

It should be remarked that the above conclusion regarding 
the validity of the “lag-in-downwash” factor differs somewhat 
from that of Martin, Diederich, and Bobbitt (ref. 28). 
These authors treat several wing-tail cases by direct considera- 
tion of the effect on the lift at the tail of the downwash shed 
by a forward wing which is plunging downward with uni- 
form acceleration aVo- For values of their results for 

Cl- and (7,„. due to interference should therefore be 

Ct ^ Ct fjp 

equivalent to those of this report for the areas — F'i{(p)d<p 

and — J F's(<p)d<p. It is found in reference 28 on compar- 
ison of some of the results with those derived from the 
simpler lag-in-downwash theory that the latter theory 
generally gives accurate results when the axis of rotation is 
near the wing centroid of area. On examination of the 
equations used for the simplified theory in reference 28, 
however, it is noted that the effective tail length (denoted 
I' in this paper) is taken to be the distance between the tail 
center of pressure and the axis of rotation, which we have 
called It- With this representation, the out-of-phase tail 
lift and moment coefficients due to interference may be 
written (in the notation of this paper) as. 




2 ^ ^ 
^ aVo 




cv 




( 83 ) 


T J 



USE OF INDICIAL FUNCTION CONCEPT IN ANALYSIS OF UNSTEADY MOTIONS OF WINGS AND WING-TAIL COMBINATIONS 37 


whereas, by the results of the present analysis, they are, 

(84) 


//O' \ o ( /O 

( M C aV, 




-2 - ^ 
int C C OcV 0 


Ct. 


The two sets of expressions agree when I' and It are equal, 
and, as noted in reference 28, the simple theory gives good 
results for just the conditions corresponding to this case. 
For other cases, however, it is believed that equations (84) 
are more representative of the physical situation, and should 
extend the range of validity of the simple theory beyond 
that indicated in reference 28. 


Another case is pointed out in reference 28 for which the 
simple theory breaks down, that is, when the forward wing 
is two-dimensional, or of large span. The results of this 
analysis tend to confirm this conclusion, since for the large- 
span wing, the down wash in steady flow is near zero. Cor- 
respondingly, the steady-state interference lift of the tail 
is near zero, which means that the relative magnitudes of the 
peaks of the detailed indicial response curve (which are 
ignored in the simple theory) play a prominent role in 

determining the magnitude of the area J Fi{(p)d<p. 


PART III 

EFFECT OF NONLINEARITIES 


In this final section we depart from discussion of the use of 
the indicial-f unction concept in dynamic-stability analyses 
to consider a problem which may have pertinence to per- 
formance characteristics under certain nonlinear flow con- 
ditions. The problem is: given a flight condition in which 
the variation of the airfoil’s static restoring moment is non- 
linear with respect to the angle of attack, what effect has 
the nonlinearity on the rotary damping-moment coefficient? 
The question is of especial pertinence to dynamic-stability 
problems of aircraft operating under separated -flow condi- 
tions, since it is well known that one effect of flow separation 
is often to cause significant departures from linearity of the 
static-restoring-moment curve. A common type of non- 
linear pitching-moment variation is that shown by the family 
of curves in figure 60, and it is this type which will be con- 
sidered here. 



Figure 60. — Family of nonlinear static pitching-moment curves. 


In order to attack the problem, we again require that the 
plunging velocity of the aircraft center of gravity be zero and, 
furthermore, assume that the nonlinear restoring moment 
is a function only of the instantaneous angle of attack. The 
characteristic differential equation governing the motion may 
then be written as, 

a)«-l-K(a)a=0 (85) 


where I is the moment of inertia, f (a, a) represents the damp 
ing coefficient as a nonlinear function of angle of attack and 
angular velocity, and k{o) is the restoring-moment coefficient 
as a nonlinear function of angle of attack. 

Next, we contend that the departures from linearity of the 
restoring moment curves k{ol) are caused primarily by flow 
separation, so that in the absence of separation the curves 
would be nearly linear. In any case, however, a given 
restoring-moment variation can be broken into two parts; 
one linear, and given, say, by the results of linearized-flow 
theory, the other, a nonlinear variation of such magnitude 
as to combine with the linear variation to give the observed 
curve. These statements are illustrated in figure 61. 



Figure 61. — Breakdown of pitching-moment curve into linear ami 
nonlinear components. 


It must be noted that the restoring-moment curves shown 
in figures 60 and 61 are presumed to be taken from the 
results of static wind-tunnel measurements, where the angle 
of attack is simply the angle between the chord-plane of 
the' wing and the free-stream direction. When the wing is 
oscillating, however, each point ^ measured from the axis of ro- 


a? 


tation experiences an additional angle of attack due to the 


angular velocity. Then if the steady-state pitching-moment 
variation presumed to be caused by flow separation (curve 
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(c) in figure 61) is approximated in the operating angle-of- 
attack range by, 

Cm{ot)—aoi—ho? 
in the unsteady case it is 

( 86 ) 

where now ^ is assumed to be the distance from the axis to 
the point at which the additional lift due to separation is 

O'? 

concentrated. Now, since for slow frequencies, ^ is much 

smaller than a, all but first-order terms in a are neglected, 
so that equation (86) becomes 

Cm=oi{a-ho?)-\-^{a—Ua^ (87) 

V 0 

The characteristic equation of motion then becomes 

+ j (a-3 6a^)]- a (a- M] =0 

( 88 ) 

where (C'm^+C'm^) and (C',„^) are the (constant) stability 
derivatives which would be present alone in the absence of 
separation. Thus it appears that even if the nonlinearity in 
the static restoring moment is of small magnitude, neverthe- 


less the possibility still remains that the damping moment 
can be profoundly affected in the event that is 

sufficiently small. Notice that when (Cm^+Cm^)+^^ is 

greater than zero, and if the nonlinear term in the restoring 
moment can be ignored in equation (88), that equation then 
takes on the form, 

a — 25a(l — ij.a^)-{-K^a = 0 (89) 

Equation (89) is then recognized as being the well-known 
Van Der Pol equation of nonlinear mechanics. It is evident 
that for small values of a the damping term is negative, lead- 
ing to a divergent oscillation, whereas for larger a the damp- 
ing term is positive. A stable regime therefoi'e will exist 

near a.=^> and oscillations of either large or small ampli- 

tude will converge to that regime. 

Hence, the result of this analysis would appear to indicate 
that nonlinearities in the variation of restoring moment with 
angle of attack of the type shown in figure 60 may have the 
effect of promoting self-sustained rotary oscillations of small 
amplitude in cases where the damping moment existing 
under unseparated flow conditions is small. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 16, 195 Jf. 



APPE^ 

RESPONSE IN LIFT OF TWO-DIMENSIONAL T 

The results for the section lift response to the two-dimen- 
sional vortex field are listed below, together with the limits 
of the intervals of “time” t within which the results are 
applicable. The intervals have been illustrated graphically 
in figure 38 of the text. 


Interval 1, 


l—<t< 

M,+l- - M„-M 
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’ M„4-l 

Hi da— r Il2da 
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APPENDIX B 


RESPONSE IN LIFT OF TWO-DIMENSIONAL, RECTANGULAR, AND WIDE TRIANGULAR TAILS TO TWO-DIMENSIONAL VORTEX 

SYSTEM— GUST ANALYSIS 


Approximate solutions for the responses of two-dimension- 
al, rectangular, and wide triangular tails are listed below. 
Solutions were obtained by means of a gust tj^pe analysis, 
explained in the text. The indicial gust functions used in 
the analysis are those obtained from references 3, 5, and 6. 

TWO-DIMENSIONAL TAIL 


Interval 1, 


^ <t< 


M„+l - - Mo+l 


<T 


Interval 2, <t< 


M„+l - - 
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xj c 
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i 


— — lo^ 
Mo L 
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XxfEIEwl, JL cos- 
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ir^VoC 


(--?)} 
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Interval 1, 


Mo+l 


RECTANGULAR TIP REGIONS 
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Mn+ 1 Mo 1 


<t 




where, 

V?hp=w]} 


and 


T = t 

n—t 


Mo + 1 

c 

"MX 

2A$ 


+VoC 

The response for the rectangular tail may be formed by 
direct addition of the solutions given above. 
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WIDE TRIANGULAR TAIL 
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APPENDIX C 


KESPONSE IN LIFT OP APEX-FORWARD AND APEX-REARWARD WIDE TRIANGULAR TAILS TO TWO-DIMENSIONAL VORTEX 

SYSTEM 


Exact solutions for the responses of apex-forward and 
apex-rearward wide triangular tails are listed below. Solu- 
tions were obtained by the application of strip theory. 
The section-lift response used in the analysis has already 
been presented as Appendix A. 
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The integral I n„da may be evaluated easily, and gives, 
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APEX-FORWARD TRIANGULAR TAIL 

The form and intervals of the results for the apex-rearward 
tail also apply to the apo.x-forward tail. The terms !!„, 
Ill, II 2 , are, of course, different, and these are listed below. 
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APPENDIX D 

BOUNDARY CONDITIONS AT THE TAIL DUE TO PENETRATION OF VELOCITY FIELD OF TRAILING-VORTEX SYSTEM 


The boundary conditions at a spanwise section of the tail 
due to the tail’s penetration of the normal-velocity field 
created by a trailing-vortex loop are illustrated below. The 
five types of conditions are those possible at various sections 
of the tail due to its shape and its position with respect to 
the moving and trailing vortices. 

Case 1 , 



Mg{l+^+c)+ -^{l+^+cY-^y 


ci{U)= — 


4 A$ 






yS 2tVi 

-1 % 


^cos ' 


? V(/+ 0^- W+ 

_§y 


Case 2, 
or 


cos- 


Mo |3 


)] 



ti=y 

, Mo(l+^+c)+^/{l+^+cr-fiY 




4 1 

/3 2irVoC y 


i+i 


[V(^ + ^ V(m 7 =W + 

)] 


/Sy^cos cos ~‘ — — 




Case 2(a), 


M„ /3 



M„(i-t-g+c)- V(i+?+c)^-/3y 


/3^ 

0 


, M„(^+^-{-c)+V^?+c)’’-/Sy 

' 2 — 


,, , 4 A# 


/3 27tF„c y 

( 


i [V(m+y=w^- V(^+ ^)^-/sy -b 

^ + ?-bc)] 


iSy(cos ^ 7 ^^— cos-’ — ^ 


Case 3, 




~r^«^sr 




USE OP INDICIAL FUNCTION CONCEPT IN ANALYSIS OP UNSTEADY MOTIONS OP WINGS AND WING-TAIL COMBINATIONS 43 


i\=y 


Moa+?+c)+V(^+i-hc)^-^V 

0^ 


C,(<2)=- 

Case 4, 
or 


I cos 

' ' ; ■ Mo 


-1 


__§y_~ 

l+^ + C_ 



Mo(l + ^ + c)-y'(l+^ + cy - <3 V 


. Mo{l + ^ + c)+y{l + k + cf-0^f 

h- 

^\c lb ri^] 

REFERENCES 

1. Zimmerman, Charles H.: An Analysis of Longitudinal Stability 

in Power-Off Flight with Charts for Use in Design. NACA 
Rep. 521, 1935. 

2. Wagner, Herbert; Uber die Entstehung des dynamischen Auf- 

triebes von Tragflugeln. z.f.a.M.M. Bd. 5, Heft 1, Feb. 1925, 
S. 17-35. 

3. Heaslet, Max A., and Lomax, Harvard; Two-Dimensional Un- 

steady Lift Problems in Supersonic Plight. NACA Rep. 945, 
1949. 

4. Lomax, Harvard, Fuller, Franklyn B., and Sluder, Loraa; Two- 

and Three-Dimensional Unsteady Lift Problems in High-Speed 
Flight. NACA Rep. 1077, 1952. (Supersedes NACA TNs 
2403 and 2387). 

5. Miles, John W.; Transient Loading of Wide Delta Airfoils at 

Supersonic Speeds. NAVORD Rep. 1235, 1950. 

6. Miles, John W.; Transient Loading of Supersonic Rectangular 

Airfoils. Jour. Aero. Sci., vol. 17, no. 10, Oct. 1950, pp. 647-652. 

7. Garrick, I. E.; On Some Reciprocal Relations in the Theory of 

Non-stationary Flows. NACA Rep. 629, 1938. 


8. Jones, Robert T.; The Unsteady Lift of a Wing of Finite Aspect 

Ratio. NACA Rep. 681, 1940. 

9. Churchill, R. V.; Modern Operational Mathematics in Engineering. 

McGraw-Hill Book Company, N. Y., 1944. 

10. Miles, John W.; The Application of Unsteady Flow Theory to the 

Calculation of Dynamic Stability Derivatives. North American 
Aviation Rep. AL-957, 1950. 

11. Jahnke, Eugen, and Emde, Fritz; Tables of Functions with For- 

mulae and Curves. Fourth ed., Dover Pub., N. Y., 1945. 

12. Durand, W. F., ed.; Aerodynamic Theory. Vol. V, Julius Springer 

(Berlin), 1935. (Available as CIT Reprint, 1943). 

13. Theodorsen, Theodore; General Theory of Aerodynamic Instabil- 

ity and the Mechanism of Flutter. NACA Rep. 496, 1935. 

14. Garrick, I. E., and Rubinow, S. I.; Flutter and Oscillating Air- 

Force Calculations for an Airfoil in a Two-Dimensional Super- 
sonic Flow. NACA Rep. 846, 1946. 

15. Watkins, Charles E.; Effect of Aspect Ratio on the Air Forces and 

Moments of Harmonically Oscillating Thin Rectangular Wings in 
Supersonic Potential Flow. NACA TN 2064, 1950. 

16. Watkins, Charles E.; Air Forces and Moments on Triangular and 

Related Wings with Subsonic Leading Edges Oscillating in 
Supersonic Potential Flow. NACA TN 2457, 1951. 

17. Nelson, Herbert C.; Lift and Moment on Oscillating Triangular 

and Related Wings with Supersonic Edges. NACA TN 2494, 
1951. 

18. Miles, John W.; On Harmonic Motion of Wide Delta Airfoils at 

Supersonic Speeds. NAVORD Rep. 1234, 1950. 

19. Ribner, Herbert S., and Malvestuto, Frank S., Jr.; Stability De- 

rivatives of Triangular Wings at Supersonic Speeds. NACA 
Rep. 908, 1948. 

20. Malvestuto, Frank S., Jr., and Margolis, Kenneth; Theoretical 

Stability Derivatives of Thin Sweptback Wings Tapered to a 
Point with Sweptback or Sweptforward Trailing Edges for a 
Limited Range of Supersonic Speeds. NACA Rep. 971, 1950. 
(Supersedes NACA TN 1761). 

21. Malvestuto, Frank S., Jr., and Hoover, Dorothy M.; Lift and 

Pitching Derivatives of Thin Sweptback Tapered Wings with 
Streamwise Tips and Subsonic Leading Edges at Supersonic 
Speeds. NACA TN 2294, 1951. 

22. Malvestuto, Frank S., Jr., and Hoover, Dorothy M.; Supersonic 

Lift and Pitching Moment of Thin Sweptback Tapered Wings 
Produced by Constant Vertical Acceleration. Subsonic Leading 
Edges and Supersonic Trailing Edges. NACA TN 2315, 1951. 

23. Jones, Robert T., and Fehlner, Leo F.; Transient Effects of the 

Wing Wake on the Horizontal Tail. NACA TN 771, 1940. 

24. Lomax, Harvard, Heaslet, Max A., and Fuller, Franklyn B.; 

Three-Dimensional, Unsteady-Lift Problems in High-Speed 
Flight-Basic Concepts. NACA TN 2256, 1950. 

25. Miles, John W.; Simple Planforms in Supersonic Flow. Jour 

Aero. Sci., vol. 17, no. 2, Feb. 1950, p. 127. 

26. Walsh, J., Zartarian, G., and Voss, H. M.: Generalized Aerody- 

namic Forces on the Delta Wing with Supersonic Leading Edges. 
Jour. Aero. Sci. vol. 21, no. 11, Nov. 1954, pp. 739-748. 

27. Heaslet, Max A., and Spreiter, John R.; Reciprocity Relations in 

Aerodynamics. NACA Rep. 1119, 1953. 

28. Martin, John C., Diederich, Margaret S., and Bobbitt, Percy J.; 

A Theoretical Investigation of the Aerodynamics of Wing-Tail 
Combinations Performing Time-Dependent Motions at Super- 
sonic Speeds. NACA TN 3072, 1954. 

29. Ribner, Herbert S.; Time-Dependent Downwash at the Tail and 

the Pitching Moment due to Normal Acceleration at Supersonic 
Speeds. NACA TN 2042, 1950. 


U. S. GOVERNMENT PRINTING OFFICE: 1955 



